AD  No.- — 

DDC  FILE  COPY;  ADA043377 


02  000173.00 


THE  THEORETICAL  BASIS 
OF  THE  CODE  50 
NUCLEAR  EXCHANGE  MODEL 

Neal  H.  Hillerman 


Research  Contribution  173 


I 


Center 
for 
NavaL 
Analyses 


Systems  Evaluation  Group 


an  tffilittt  of  tht  UniYtnity  of  Rochtttar 
1401  Wilton  Boubvard,  Arlington,  Virginia  22309 


NOTICE 


THIS  DOCUMENT  YA3  BS£N 
AWSftW  FOR,  PUBLIC  RELEASE, 

. - - »*  »|  '/.  ‘.VA.VhK'Uit 


^ AUG  ^ 

teraranns/yi 


>(  ^ 


/■  r o • 
a.,  . 


■'V' 


T 


1401  Wilson  Boulevard  Arlington.  Virginia  22200 


^ A A A|<1.[|A| 


703/524-9400 


Center 
for 

-**aa**^/  dl^d/ 


Analyses 


(SEG) 304-72 
21  Mar  1972 


From:  Director,  Systems  Evaluation  Group 

To:  Distribution  List 


Subj:  Center  for  Naval  Analyses  Research  Contribution  173- 

forwarding  of 


Enel:  (1)  CRC  173,  "The  Theoretical  Basis  of  the  code  50 

Nuclear  Exchange  Model  (U) , " by  N.H.  Hillerman 
September  1971,  Unclassified 


1.  Enclosure  (1)  is  forwarded  as  a matter  of  possible  interest. 

2.  Although  the  code  50  model  described  in  this  research 

tS  n0t  thS  m°St  uP~to“date  version  available,  it 
**  P blt  hS<?  provide  documentation  of  the  theoretical  and 
mathematical  background  upon  which  all  code  50  models  are  based 


valu^in^  COn^rJbutions  are  distributed  for  their  potential 

reorLent  ^er  S and  analyses-  They  do  not  necessarily 

represent  the  opinion  of  the  Department  of  the  Navy. 


,/Wi  e- ' ^ 

PHIL  E.  Deppy 
Director 


Distribution  List: 
Reverse  page 


An  Equal  Opportunity  Employer 


(SEG) 304-72 
21  Mar  1972 


Subj : Center  for  Naval  Analyses  Research  Contribution  173? 

forwarding  of 


DISTRIBUTION  LIST: 

Department  of  Defense 
JCS  (SAGA) 

DDR&E  DIR,  DNA 

NAVMBR,  STP  COMDT,  AFSC 

SENAVMBR,  WSEG  COMDT,  ICAF 

Department  of  State: 

U.S.  Arms  Control  and  Disarmament  Agency 

Department  of  the  Navy: 

SECNAV 

OPA 

ONR 

OpNav : Op-09B9,  Op-96,  Op-06,  Op-604,  Op-963G 

CINCPACFLT 

Department  of  the  Air  Force: 

Dept  of  AF  (Attn:  PDO  4008) (9) 
Superintendent  Air  University 
Air  Proving  Ground  Center,  Eglin  AFB 

Other: 

ANSER 
APL.  .1U 
IDA 
RAND 
SRI 

Cornell  Aeronautical  Lab,  Inc. 

Sandia  Corp. 

Aerospace  corp. 

Lambda  corp. 


CNM  PM-1,  Strategic  Systems  Project  Office 


CINCLANTFLT 

CINCUSNAVEUR 

NUWPNTRAGRULANT 

NUWPNTRAGRUPAC 

COMOPTEVFOR 

DEPCOMOPTEVFORPAC 

CINCPAC 

CINCLANT 

CINCSAC 

DIR,  NRL 

STRATANALSUPPGRU 

DODCI 

FLECOMPUTPROGCENLANT 

FLECOMPUTPROGCENPAC 

USNA  ANNA 

NAVPGSCOL 

NAVWARCOL 

NAVAI RDEVCEN 

NAVWPNSCEN 

NWL 

NOL 

NAVWPNSYSANALO 
NAVAI RSYSCOMHQ 
NAVORDSYSCOMHQ 


CENTER  FOR  NAVAL  ANALYSES 


Systems  Evaluation  Group  ^ 

Research  Contribtrttow-173  j 


THE  THEORETICAL  BASIS  OF  THE 
CODE  50'NUCLEAR  EXCHANGE  MODEL, 


CLEARED 

FOR  OPEN  PUBLICATION. 

JUL  1 3 1977  11 

JIRECTORAIE  run  ■ titrURMATIOfl 

AMO  SECURITY  REVIEW  (OASD-PA) 

O'-'’ 


This  Research  Contribution  does  not  necessarily  represent 
the  opinion  of  the  Department  of  the  Navy. 


8-A-^oei if 


Work  conducted  under  contrac^N0.0pi4— 08— A— ^0091 

y 


.Enclosure-f  1 j'tcTfSEG  l304-72dated  2 1 Marchl972 


/ , 


^/(J  . ^ '>C 


IC  .*  . '.  cUon  rtf 


J U SeUioll  □ 

□ 


Pv 


EijiK’AlI’J.Vli'f.K'Il'f  M I 

. SPECIAL 


1 


y 


L 


ABSTRACT 

"^The  Code  50  Nuclear  Exchange  Model  is  a war  game  model  produced  by  the  LAMBDA 
Corporation.  This  Research /Contribution  derives  and  explains  the  basic  mathematical 
models  used  in  the  computer  programs  of  that  model,  including  models  of  missile  and 
bomber  penetration  as  well  as  damage  calculation,  weapon  allocation,  and  kill  probability 
models.  Model  implementation  and  integration  into  the  Code  50  program  are  also 
demonstrated. 
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SECTION  I 


INTRODUCTION 


GENERAL 

This  Research  Contribution  presents  the  derivations  and  explanations  of  the  basic 
mathematical  models  used  in  the  computer  programs  of  the  Code  50  nuclear  exchange 
model  and  includes  modifications  by  CNA  to  permit  exercise  of  the  Code  50  model  and 
submodels  on  computers  of  limited  core -storage  capacity.  The  models  included  in  the 
Code  50  nuclear  exchange  model  are  concerned  with: 

• Missile  penetration 

• Bomber  penetration 

• Damage  calculation 

• Weapon-target  allocation 

• Single -shot  kill  probability. 

The  final  section  of  this  Research  Contribution  shows  how  the  various  models  are 
implemented  and  integrated  into  the  Code  50  program.  Additional  documentation  on  the 
exact  details  of  computer  operations  of  the  Code  50  program  is  published  in  Research 
Contribution  133,  The  Calculation  Procedure  for  the  CNA  Version  of  the  Code  50  Nuclear 
Exchange  Model . '">ad  37^ 

NUCLEAR  EXCHANGE  MODELS 

A nuclear  exchange  model  simulates  the  result  of  strategic  warfare  in  a ?. -player 
game  in  which  nuclear  weapons  only  are  employed.  There  are  basically  two  types  of 
such  models: 

The  first  includes  highly  detailed  simulations,  used  to  establish  targeting  strategy 
for  various  contingencies.  Such  simulations  consider  many  variables,  requiring,  in  some 
cases,  weeks  of  data  preparation  and  many  hours  of  computer  time. 

The  second  type  includes  the  many  aggregated  gross  effectiveness  models  in  use 
today.  In  this  type  of  model,  many  individual  targets  are  grouped  in  a single  target  class 
if  their  individual  characteristics  are  similar.  Such  details  as  range  constraints,  attack 
azimuth,  attack  timing,  etc.,  are  not  considered.  The  Code  50  exchange  model  is  of  the 
latter  type. 

In  the  play  of  a nuclear  exchange  war  game,  it  is  first  necessary  to  decide  which  side 
strikes  first,  second,  third,  etc.,  and  how  many  strikes  each  delivers.  It  is  also  neces- 
sary to  specify  whifch  types  of  target  are  attacked  in  each  strike.  For  example,  targets 
might  be  the  strategic  force,  population  centers,  or  mixes  of  n .itary  targets  and  popula- 
tion centers.  The  terms  "counter -force;**  "counter-value,"  and  "mixed"  are  used, 
respectively,  to  define  these  types  of  target. 
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In  a counter-force  attack,  the  attacker  attempts  to  minimize  future  damage  to  him- 
self by  striking  at  the  other  side’s  strategic  nuclear  force.  This  is  also  known  as  a 
" damage  -limiting”  attack . 

In  a counter-value  strike,  the  attacker  attempts  to  maximize  the  damage  to  some 
measure  of  value  on  the  other  side.  Generally,  this  measure  of  value  is  expressed  in 
terms  of  either  population  or  industrial  floor  space.  Mixed  targets  present  some  of  the 
characteristics  of  both  counter-value  and  counter-force  targets. 

Since  mixed  targets  present  complications,  particularly  with  respect  to  target  value 
assignments,  the  Code  50  nuclear  exchange  model  considers  pure  counter-force  and 
counter-value  strikes  only. 

CODE  50  BACKGROUND 

Code  50  is  an  early  version  of  a nuclear  war  game  developed  by  the  Lambda  Corpora- 
tion for  the  Office  of  the  Assistant  Secretary  of  Defense  (Systems  Analysis).  Throughout 
this  model,  expected  values  that  tend  to  be  conservative  from  the  point  of  view  of  the 
offense  ("offense -conservative")  are  used.  For  example,  the  probability  of  penetrating 
defenses  tends  to  be  offense -conservative.  Blast  damage  to  the  intended  target  is  the 
only  damage  considered,  with  the  result  that  both  collateral  damage  and  damage  resulting 
from  effects  other  than  blast  are  not  included.  Hence,  total  damage  inflicted  also  tends 
to  be  offense -conservative. 

Code  50  allows  a maximum  of  three  strikes . The  choice  of  the  first  attacker  is 
determined  arbitrarily;  on  the  second  strike,  however,  the  roles  of  the  attacker  and 
defender  are  reversed.  The  roles  are  reversed  again  on  the  third  strike.  The  type  of 
strike  is  also  arbitrary  on  the  first  two  strikes  (i.e.,  counter-value  or  counter-force), 
but  the  third  strike  is  always  a counter-value  strike. 

The  forces  and  cities  of  the  two  sides  are  agglomerated,  each  side  being  restricted 
to  a total  of  48  target  types . These  include  all  the  weapon  types  and  all  the  city  types  for 
the  side.  The  program  does  not  require  that  the  proportion  of  weapon  types  to  city  types 
conform  to  any  pre-established  ratio. 

CODE  50  EXPERIENCE  AT  CNA 

Under  Task  1 of  the  Strategic  Force  Study,  the  Center  for  Naval  Analyses  was 
charged  with  evaluating  Code  50.  This  included  program  adaptation  to  exercise  the 
Code  50  model  on  the  CNA  computer,  as  well  as  documentation  and  evaluation  of  the 
program.  Because  the  original  Code  50  program  required  more  core  storage  than  the 
CNA  computer  could  provide,  modifications  were  made  to  reduce  the  core  requirements. 
The  resulting  program,  designated  "CNA  50",  carried  out  generally  the  same  operations 
and  yielded  approximately  the  same  results  as  earlier  implementations  of  the  original 
Code  50  program. 
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SECTION  II 


CODE  50  MATHEMATICAL  MODELS 


GENERAL 

This  section  explains  and  derives  the  appropriate  mathematical  expressions  for 
each  of  the  Code  50  models,  to  promote  a fuller  understanding  of  the  underlying  implica- 
tions of  the  models  themselves  and  to  provide  greater  insight  into  the  meaning  of  the 
overall  results  obtained  when  Code  50  models  are  used. 

In  both  the  missile  and  bomber  penetration  models,  the  derivation  is  in  two  parts. 
The  first  is  concerned  with  penetration  of  area  defenses;  the  second  is  concerned  with 
the  problem  of  penetrating  terminal  defenses.  In  the  case  of  bombers,  the  area  defenses 
consist  of  waves  of  fighters  sent  out  to  intercept  and  shoot  down  the  bombers.  In  the  case 
of  missiles,  the  area  defenses  consist  of  a nationwide  area  ABM  defense.  Terminal  de- 
fenses against  bombers  consist  of  surface-to-air  missiles  (SAMs);  terminal  defense 
against  missiles  is  assumed  to  be  provided  by  a terminal  ABM  system  around  the  target. 
An  additional  complication  is  provided  in  the  case  of  bombers:  those  which  penetrate 
terminal  defenses  are  allowed,  in  certain  cases,  to  attack  the  SAM  sites  directly,  rather 
than  merely  suffer  attrition  in  exhausting  the  supply  of  SAMs. 

The  single-shot  kill  probability  model  covers  the  mechanics  of  kill,  taking  into  ac- 
count target  hardness  and  weapon  yield.  The  model  considers  the  effects  of  probabilistic 
delivery  errors  and  the  geographic  distribution  of  an  individual  target’s  value,  under  the 
assumption  that  both  delivery  error  and  target  value  can  be  approximated  by  circular 
normal  distributions . 

The  damage  assessment  models  recognize  two  types  of  targets: 

The  first  is  a point  target,  that  is,  a target  with  no  appreciable  area.  Weapons  and 
other  military  targets  are  considered  point  targets  within  the  meaning  of  the  program. 
Damage  to  point  targets  is  assessed,  using  an  exponential  law  (Bernoulli  trial  methodol- 
ogy), since  each  weapon  is  assumed  capable  of  destroying  the  target  if  placed  correctly. 

The  second  type  is  an  area  target.  All  cities  are  treated  as  area  targets.  For  area 
targets,  the  square  root  law  is  employed,  since  it  is  assumed  that  a single  weapon  may 
not  be  able  to  destroy  the  entire  target,  even  if  placed  optimally. 

The  weapon -target  allocation  model  in  Code  50  has  a two -fold  purpose:  to  inflict 
the  maximum  amount  of  damage  on  the  whole  target  base  and  to  inflict  the  required 
amount  of  damage  on  preferred  targets.  Preferred  targets  are  defined  as  non -military 
targets  that  have  terminal  defenses.  The  generalized  LaGrange  multiplier  technique  is 
used,  along  with  some  heuristic  non-linear  programming  logarithms  to  improve  the 
efficiency  of  the  program  in  finding  a new  optimum  solution  that  satisfies  the  maximiza- 
tion problem  and  the  preferred  target  damage  constraint. 

MISSILE  PENETRATION  MODELS 

Offensive  missile  penetration  in  Code  50  consists  of  a two-part  problem.  First, 
there  is  the  problem  of  missiles  penetrating  an  area  defense  system  that  is  assumed 
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capable  of  protecting  (up  to  the  level  of  available  interceptors)  all  probable  targets . The 
second  problem  is  the  penetration  of  terminal  defenses  at  the  individual  targets  by  mis- 
siles that  succeed  in  penetrating  the  area  defenses . In  this  subsection,  the  area  and 
terminal  penetration  models  used  in  Code  50  are  described,  beginning  with  the  basic 
area  penetration  model  and  proceeding  to  the  terminal  penetration  and  entry  price  models. 
Subsequently,  sample  calculations  are  made  to  demonstrate  the  operation  of  the  models. 

The  subscript  conventions  in  this  subsection  are  as  follows: 

IW  Weapon  type  subscript  (missile  types  only) 

KS  Strike  subscript 

To  simplify  notation,  no  target  type  subscripts  are  used  in  this  development.  However, 
because  the  equations  do  not  vary  with  target  type,  no  confusion  should  result. 

Other  symbols  used  in  the  equations  are  listed  below: 

AIN  Number  of  interceptors  used  on  a given  strike 

AINX  Total  number  of  area  interceptors  available 

AMP™  Price  in  missiles  successfully  launched  which  just  exhausts  area 
interceptors 

AOBjw  Number  of  area  objects  per  type  IW  missile 

ARR™  „„  Number  of  missiles  of  type  IW  successfully  launched  and  arriving  at 
’ defender's  area  defenses  on  strike  KS 

ARR*  Number  of  missiles  arriving  at  area  defenses,  just  exhausting  area 

defenses  but  resulting  in  no  payoff 

ARRMX™  Maximum  number  of  type  IW  missiles  that  could  be  successfully 

launched  by  attacker  and  arrive  at  defender's  area  defenses  on  first 
strike 

FSENTjyy  ^Fraction  of  missiles  of  type  IW  sent  on  strike  KS 

NLEFTjyy  Number  of  offensive  missiles  of  type  IW  from  original  inventory 
that  have  been  withheld  from  prior  strikes 

NMINT  Number  of  terminal  interceptors 

NNWP^  Number  of  offensive  missiles  of  type  IW  in  attacker’s  inventory 

NPj^  Price  in  terms  of  independent  targetable  weapons  successfully  launched 

NWP^y  Total  number  of  independent  targetable  weapons  available 

OARR^y  Number  of  area  objects  successfully  launched  by  the  attacker  and 

w arriving  at  defender's  area  defenses  on  strike  KS 

OBJ™  Maximum  number  of  area  objects  that  could  be  launched  by  the  attacker 

l and  arrive  at  the  defender’s  area  defenses  on  the  first  strike 
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OBJD 

OPEN 

PENIW 

PENPROB 

PKMIS 


RELAI 

RNIW,KS 

rriw,ks 

suRVjyv 


TOB 

TOS. 


IW 


'IW 
TPRWPN 


IW 


Number  of  area  objects  that  can  be  destroyed  by  area  interceptors 
Number  of  area  objects  penetrating  the  area  defense 
Equivalent  number  of  type  IW  missiles  penetrating  area  defense 

Probability  that  missile  area  objects  will  penetrate  the  area  defense 

Probability  of  kill  of  an  individual  interceptor  against  an  individual 
area  object 

Interceptor  reliability 

Probability  that  missile  type  IW  does  not  fail  during  flight  on  strike  KS 

Probability  that  missile  type  IW  does  not  fail  during  launch  on  strike  KS 

Probability  of  NLEFT^y  missiles  surviving  prior  strikes  by  enemy 

Terminal  objects  per  missile 
Number  of  terminal  objects  surviving 

Number  of  independently  targetable  weapons  per  missile  for  weapon 
type  IW 


Assumptions  Concerning  Penetration  of  Area  Defenses 

Basic  assumptions  regarding  area  defenses  are  as  follows: 

• Area  defenses  are  capable  of  defending  all  targets  and  may  not  be  bypassed. 

• Area  defenses  are  so  balanced  that  there  is  no  more  advantage  in  attacking 
the  area  defense  weapons  than  in  attacking  targets  directly. 

• The  area  defenses  do  not  preferentially  defend  but  randomly  destroy  the 
incoming  offensive  weapons . 

• The  defense  is  unable  to  discriminate  between  warheads  and  decoys. 

• Offensive  weapon  spacing  is  assumed  such  that  a single  interceptor  may 
destroy  only  a single  area  object. 

• On  the  first  strike,  the  defender's  strategy  is  to  use  only  a fraction  of  his 
area  interceptors;  that  fraction  is  equal  to  the  ratio  of  area  objects  actually  arriving  to 
the  number  of  area  objects  that  could  arrive. 

• The  offense  knows  exactly  how  many  of  its  weapons  have  survived  prior 

strikes. 

• Losses  because  of  unreliability  occur  before  the  missile  area  objects  reach 
the  area  defenses  and  do  not  draw  down  the  area  interceptors. 

The  Area  Defense  Penetration  Model 

The  significant  parameter  in  the  area  defense  penetration  model  is  PENPROB,  the 
probability  that  an  individual  independent  weapon  in  the  remaining  inventory  will  pene- 
trate the  area  defense. 
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Initially,  there  are  NNWPp^  offensive  missiles  of  type  IW  in  the  attacker's  inven- 
tory. At  present,  NLEFT^y  offensive  missiles  of  type  IW  from  the  original  inventory 
have  been  withheld  from  prior  strikes.  The  probability  that  these  NLEFT^  missiles 
survived  prior  strikes  by  the  other  side  is  SURV^  . The  conditions  on  the  three  strikes 
are: 


NLEFT^  = NNWP^ 

SURVjyy  =1 
NLEFTIW  = NNWPIW 
SURV^  si 

NLEFT^  = NNWP^  (1-  FSENT^  ^ 
SURV1W  s 1 


KS  = 1 


KS  = 2 


KS  = 3 


U) 


The  reprogrammable  reliability  (the  launch  reliability  or  the  probability  of  failing 
such  that  remaining  weapons  can  be  retargeted,  taking  the  failure  into  account)  for 
missile  type  IW  on  strike  KS  is  RR^  £g  . The  non -reprogrammable  reliability 

(in-flight  reliability  or  the  probability  of  not  failing  after  it  is  too  late  to  retarget,  taking 
the  failure  into  account)  for  missile  type  IW  on  strike  KS  is  RN^  Kg  . The  fraction 

of  the  remaining  missiles  of  type  IW  sent  on  strike  KS  is  FSENT^  Kg  . 


The  maximum  number  of  missiles  of  type  IW  that  could  be  successfully  launched 
by  the  attacker  and  arrive  at  the  defender’s  area  defenses  on  the  first  strike  is: 

ARRMX^  = NNWP^  • RR^  l • RN^  l (2) 


Now,  for  missile  type  IW  there  are  AOB^y  area  objects  per  missile.  Therefore, 

the  maximum  number  of  area  objects  that  could  be  launched  by  the  attacker  and  arrive 
at  the  defender's  area  defenses  on  the  first  strike  is: 

OB]  =Z  ARRMX^  * AOB^  (3) 

The  number  of  missiles  of  type  IW  successfully  laurched  and  arriving  at  the  defen- 
der's area  defenses  on  strike  KS  is: 

ARRIW  = NLEFT^  . SURV^  ♦ FSENTj^  ^ . RR^^g  • RN^^g  (4) 

The  number  of  area  objects  successfully  launched  by  the  attacker  and  arriving  at  the 
defender's  area  defenses  on  strike  KS  is  then: 
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(5) 


OARR  =EARRw  • AOBm 

On  the  first  strike,  the  area  interceptors  are  "shared  out,"  so  that  the  fraction  of 
area  interceptors  expended  is  identical  to  the  ratio  of  offensive  area  objects  arriving 
in  strike  one  to  the  maximum  number  of  offensive  area  objects  that  could  arrive.  Thus, 
for  any  remaining  interceptors,  AIN  will  be  expended  on  strike  one  according  to  the 
relationship: 


AIN  _ OARR  ^ , 

XTNX DBF  : " 1 


(6) 


The  number  of  area  interceptors  expended  on  a given  strike  is,  then: 

AIN  = AINX  • OARR/OBJ  : KS  = 1 (7) 

AIN  = AINX  : KS  > 1 

The  number  of  area  interceptors  remaining  is,  then: 

AINX  -4 AINX  - AIN  (8) 

The  destruction  of  area  objects  by  area  interceptors  is  determined  on  the  basis  of 
zero  leakage  until  the  number  of  area  objects  exceeds  the  expected  kill  potential  of 
the  area  interceptors.  Thereafter,  the  numbei  of  area  objects  destroyed  is  random 
and  proportional  to  the  number  of  interceptors  (AIN),  interceptor  reliability  (RELAI), 
and  the  probability  of  kill  of  an  individual  interceptor  against  an  individual  area  object 
(PKMIS).  The  maximum  number  of  area  objects  that  can  be  destroyed  is  therefore: 


OBJD  = AIN  • RELAI  • PKMIS 


(9) 


The  number  of  area  objects  penetrating  the  area  defense  is,  then: 

OPEN  = MAX  < 0 : OARR  - OBJD  > (10) 


The  probability  that  missile  area  objects  penetrate  the  area  defense  is,  then: 


PENPROB= 


OPEN 

iMir 


(ii) 


The  equivalent  number  of  missiles  of  a given  type  penetrating  the  area  defense,  then, 
3s: 


PEN^  = ARR^  • PENPROB  (12) 

Terminal  Defense  Penetration  and  Total  Price  to  Attack  a Target 
These  are  the  assumptions  regarding  terminal  defenses: 
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A.  The  number  of  terminal  intercepts  is  equal  to  the  total  number  of  terminal 
interceptors.  (This  implies  either  that  each  terminal  interceptor  is  perfectly  reliable  or 
that  a reduced  number  of  interceptors  is  assigned  to  compensate  for  reliability  losses.) 

B.  Terminal  interceptors  have  a probability  of  kill  of  unity  against  all  terminal 
objects  arriving  at  the  target. 

C.  The  terminal  defense  may  not  be  suppressed,  and  no  terminal  objects  penetrate 
to  the  target  until  the  terminal  interceptor  supply  is  exhausted. 

D.  The  defense  is  unable  to  discriminate  between  warheads  and  decoys. 

E.  Offensive  weapon  spacing  is  assumed  such  that  a single  interceptor  destroys 
exactly  one  terminal  object. 

F.  All  offensive  terminal  objects  liave  perfect  reliability,  all  unreliability  losses 
having  occurred  before  penetration  of  area  defenses. 

As  presented  earlier,  the  number  of  missiles  that  penetrate  the  area  defenses  .s: 

PEN^  = ARR^  • PENPROB  (13) 

Based  on  the  assumption  of  TOB^  terminal  objects  per  missile,  NMINT  terminal 

interceptors  with  unity  kill  probability  and  perfect  reliability  can  destroy  NMINT 
offensive  terminal  objects . 

The  number  of  terminal  objects  surviving  is,  then: 

TOS^  = MAX  <0  : (ARR^  • PENPROB  . TOB^)  - NMINT  > (14) 

Now,  if  TOSjyy  is  exactly  zero,  the  number  of  missiles  arriving  at  the  area  defenses 
that  just  exhausts  the  area  interceptors  but  results  in  no  payoff  is  ARR^,*  . 

ARR^*  = 'NMINT/ (PENPROB  • TOBRV)  (15) 

Another  way  of  expressing  the  price  is  in  terms  of  missiles  successfully  launched. 

If  equation  (15)  is  divided  by  RNj^  ^g  , the  price  will  be  in  these  terms: 

AMP1W  = ARRIW*/RNIW,  KS 

Finally,  if  this  result  is  multiplied  by  the  number  of  independently  targetable  weapons 
per  missile,  and  the  result  is  converted  to  an  integer,  the  price  will  be,  in  terms  of  in- 
dependently targetable  weapons  successfully  launched: 

NPIW  = !u  CAMPIW  • TPRWPNIW  1 <i7> 


-8- 


The  total  number  of  independently  targetable  weapons  available,  N\VPp^,  is  just  the 

integer  number  of  missiles  that  will  be  successfully  launched  on  this  strike  (the  number 
left  multiplied  by  the  fraction  sent  and  the  launch  reliability),  multiplied  by  the  number 
of  independently  targetable  weapons  per  missile,  with  the  results  truncated  to  an  integer 
value. 


NVVPlW  = b 


ID[NLEFTrw  . surViw  -fsent^s  .RVks3 


. TPRWPN 


IW 


(18) 


Now,  if  NPj^,  s NWP^y,  the  price  of  entry  exceeds  or  equals  the  available  resources 
and  no  payoff  from  the  target  can  be  obtained  by  weapon  type  IW  alone. 

If  NP^y  < NWPpyy  the  entry  price  allows  a payoff  from  weapon  type  IW  . 

It  should  be  noted  that  the  preceding  calculations  of  entry  price  presumes  that  only 
a single  weapon  type  will  be  used  on  a terminally  defended  target,  since  there  are  no 
provisions  for  allocating  the  total  price  among  several  weapon  types. 

Sample  Problem  of  Missile  Penetration 

A series  of  cases  will  now  be  considered  and  the  necessary  calculationo  carried  out 
in  a tabular  manner,  with  one  case  explained  in  detail.  The  parameter  values  for  the 
cases  are  as  follows: 

KS  > 1 

NLEFT^  . SURVRV  . FSENT^  ^g  = 10,20,30,40,50 

rriw,ks  " 0-9 

RNIW,KS=  °'8 
AOBjw  = 3 

TOB^  = 3 

TPRWPN^  = 3 

AIN  = 30 

RELAI  = 0.8 

PKMIS  = 1.0 

NMINT  = 30 

The  case  of  NLEFT^y  • SURV^  • FSENT^  = 50  will  now  be  carried  out  in 
detail . 
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The  number  of  missiles  arriving  at  the  area  defenses  is: 

ARRiw  = (NLEFTiw  f SURViw  * fsentiw,ks)  * RRIW,KS  * RNIW,KS 
= (50)  • (.9)  • (.8)  = 36 

The  14  missiles  lost  consist  of  5 launch  failures  and  9 inflight  failures.  Because 
only  a single  missile  type  is  being  analyzed,  the  number  of  area  objects  arriving  is: 

OARR  = ARRj^y  ♦ AOBj^y 
= (36)  • (3)  = 108 

The  number  of  effective  area  interceptors  is: 

OBJD  = AIN  • RE  LAI  • PKMIS 
= (30)*  (.8)  • (1.0)  = 24 

Therefore,  the  number  of  area  objects  penetrating  the  area  defenses  is: 

OPEN  = MAX  <0  : OARR  - OBJD  > 

= MAX  <0  : 108  - 24>  84 


The  penetration  probability  may  then  be  found: 


PENPROB  = 


OPEN 


84 

W 


= .7778 


The  number  of  missiles  penetrating  is,  then: 


PEN^  = ARRjw  • PENPROB 
= (36)  • (.7778)  = 28 


The  number  of  terminal  objects  surviving  terminal  defenses  is: 
TOS^  = MAX  <0  : (aRR^  • PENPROB  • TOB^)  - NMINT > 
= MAX  < 0 : ((36)  • (.7778)  • (3) ) - 30  > 

= 54 


The  price  in  missiles  arriving  at  area  defenses  that  was  necessary  to  gain  a payoff 

was: 

ARRjyy*  = NMINT/(PENPROB  • TOB^) 

= 30/(. 7778  • 3)  = 12.86 
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This  translates  into  the  price  of  missiles  successfully  launched: 


AMPIW  ARR^  */RNIW  ^ KS 
= 12. 86/. 8 = 16.08 

In  terms  of  an  integer  number  of  independently  targetable  weapons,  the  price  is: 
NPIW  = IU  CAMP^  • TPRWPN^] 

= I [16.08  f 3] 

= 49 


Now,  the  number  of  allocatable  weapons  was: 

Id[NlEFtiw  • SURVlw  • FSENTiw>ks  • RRM|KS]  • TPRWPNiw 


NWPjw  - ID 


= Ir 


D 
= 135 


Id[(50)  • (.9)]  • (3) 


Since  NP^y  s NWPj^  , a payoff  may  be  secured  from  the  target. 

The  number  of  weapons  providing  the  payoff  is: 

= MAX  <0  : (NWPjW  - NP^)  • PENPROB  • RN^^g  > 

= MAX  <0  : (135  - 49)  • (.7778)  • (.8)  > 

= 54 

Notice  that  this  result  is  the  same  as  the  result  obtained  for  the  number  of  terminal 
objects  surviving  terminal  defenses  since  the  terminal  objects  per  missile  and  the 
independently  targetable  weapons  per  missile  are  both  set  to  the  same  value  (3).  The 
results  for  the  other  cases  are  shown  in  table  1. 


TABLE  I 
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SUMMARY  OF  RESULTS  FOR  SAMPLE  CASES 


Number  of 
missiles 
launched 

Number  of 
missiles 
successfully 
launched 

Number  of 
allocatable 
weapons, 
NWP,W 

Number  of 
missiles  that 
reach  area 
defense, 

ARR,W 

Probability  of  a 
missile  or  weapon 
penetrating  area 
defense, 
PENPROB 

Number  of 
missiles  that 
reach  terminal 
defense, 
PEN|W 

Number  of 
missiles  that 
arrive  at 

target 

Number  of 
targe  table 
weapons 
on  target 

10 

9 

27 

7.2 

0 

0 

0 

0 

20 

18 

54 

14.4 

.4444 

6.4 

0 

0 

30 

27 

81 

21.6 

.6296 

13.6 

3.6 

11 

40 

36 

108 

28.8 

.7222 

20.8 

10.8 

32 

50 

45 

135 

36.0 

.7778 

28.0 

18.0 

54 

BOMBER  PENETRATION  MODELS 

As  with  missile  penetration,  bomber  penetration  in  Code  50  consists  of  a two-part 
problem.  First,  there  is  the  problem  of  the  bombers  that  penetrate  the  interceptor 
defense  system.*  The  second  problem  is  penetration  of  the  SAM  defenses,**  The 
remainder  of  this  subsection  discusses  the  penetration  models,  beginning  with  the  model 
of  interceptor  defenses,  and  then  proceeds  to  the  SAM  defense  penetration  model.  Some 
sample  calculations,  showing  how  the  models  function  in  Code  50,  are  then  presented. 

The  subscript  conventions  used  in  this  subsection  are  as  follows: 

IW  Weapon  type  subscript  (bomber  types  only) 

IB  Bomber  class  subscript 

IF  Fighter  type  subscript 

KS  Strike  subscript 


*This  corresponds  to  the  penetration  of  area  defenses  in  the  missile  penetration 
problem  although,  mathematically,  the  missile  penetration  and  bomber  penetration 
are  handled  much  differently. 

**This  corresponds  to  penetration  of  the  terminal  defenses  in  the  missile  penetration 
problem . 


-12- 


«1S 


To  simplify  notation,  no  target  type  subscripts  are  used  in  this  development.  However, 
because  the  equations  derived  do  not  vary  with  target  type,  no  confusion  should  result. 


Other  symbols  used  in  the  equations  are  listed  below: 


AMP^  ’ Price  in  bombers  successfully  launched 

ARELj^,  ASM  reliability 


ARRm  Number  of  bombers  and  decoys  successfully  launched 

and  arriving  at  defender's  fighter  defense  system  on 
strike  KS 

ARRjyy*  Price  in  bombers  that  arrive  at  defenses 

ASMDECjyy  Number  of  ASM  decoys  per  ASM  warhead 


B 

bdesm 

BDPENjyy 

BMRDEC 


IW 


Total  number  of  bombers  and  decoys  arriving  at  fighter 
defense  system 

Number  of  bombers  destroyed  by  all  NSAMX  missiles 

Number  of  bombers  and  decoys  of  weapon  type  IW  expected 
to  penetrate  fighter  defenses 

Number  of  bomber  decoys  per  bomber  of  type  IW 


F 

FKBD^ 

FSENT 


IW,  KS 


IWP 


IW 


NI 


IF 
NLEFT 


IW 


Fraction  of  bombers  that  attack  through  given  corridor;  also 
fraction  of  fighters  that  attack  bombers 

Fraction  of  type  IW  bombers  killed  before  weapon  delivery 

Fraction  of  remaining  type  IW  bombers  sent  on  strike  KS 

Number  of  warheads  required  to  produce  at  least  95  percent 
kill  probability 

Number  of  fighters  in  defense  system  of  type  IW 
Number  of  bombers  that  have  not  been  sent  in  prior  strikes 


NNWP^  Number  of  type  IW  bombers  in  attacker’s  inventory 


NPSNIW 

NPSS^y 

NSAMDIT 


Total  price  in  bombers,  or  in  subsonic  cruise  armed  decoys 
(SCADs) 

Number  of  independently  targetable  weapons  successfully 
launched 

Number  of  independently  targetable  weapons  for  suppressing 
SAM  sites 

Number  of  suppressible  SAM  sites  at  target  IT 
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NSAMDM 

NSAMX 

NWP^ 

PENIW 

PENPROB 


PK 


IW 

IF'miw 

pkillt 


"IB 


PSEFF 

PSURV 

RIF 
RN 


IW 

IW 

®iw»  if 


IW,KS 

RRIW,KS 

SSK 

SSKK 

SURVm 

TPRWPN 


IW 


WHDS 


IW 


Number  of  suppressible  SAMs  per  site 

Number  of  non -suppressible  SAMs 

Number  of  independently  targetable  weapons  available 

Expected  number  of  bombers  penetrating  defenses 

Overall  penetration  probability  for  type  IW  weapon 

Probability  that  a type  IF  attacking  fighter  kills  a type  IB^ 
bomber 

Probability  that  a bomber  or  decoy  is  killed 

Probability  of  bomber  survival  until  weapon  launch 

Probability  that  none  of  the  F • NI.p  fighters  kills  a given 
bomber  of  same  type 

Ratio  of  type  IF  interceptors  to  total  number  of  bombers  in 
the  corridor 

Probability  that  bomber  type  IW  does  not  fail  during  flight 
on  strike  KS 

Probability  thar  bomber  type  IW  does  not  fail  during  launch 
on  strike  KS 

Single -shot  kill  probability  for  ASM  bomber,  taking  into 
account  SAM  site  hardness,  SRAM  yield,  CEP,  and  height 
of  burst 

Overall  kill  probibility  of  ASM  bomber  against  suppressible 
SAM  site 

Probability  that  NLEFT-,  bombers  survive  prior  enemy 
strikes 

Number  of  independently  targetable  weapons  per  missile 
for  weapon  type  IW 

Number  of  ASM  warheads  per  bomber  for  weapon  type  IW 


Interceptor  Penetration  Model 

Before  a discussion  of  the  penetration  equations  themselves,  we  shall  list  the  basic 
assumptions  regarding  the  interceptor  defenses: 

A.  The  interceptors  are  assumed  capable  of  defending  all  targets  and  may  not  be 
bypassed. 

B.  There  is  no  way  of  destroying  the  interceptors  on  the  ground  before  the  bombers 
attempt  penetration. 

C.  The  interceptors  attack  the  bomber  force  randomly  in  waves,  each  wave  con- 
sisting of  a single  type  of  interceptor. 
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D.  The  interceptor/bomber  ratio  is  constant  for  a single  interceptor  type  in  all 
attack  corridors. 

E.  There  is  no  coordination  between  interceptor  waves;  i.e.,  the  whole  force  of 
interceptors  of  each  type  attacks  the  whole  bomber  force  on  a random  basis. 

F.  The  defense  is  unable  to  discriminate  between  bombers  and  bomber  decoys 
and  between  ASMs  and  ASM  decoys . 

G.  The  offense  knows  exactly  how  many  of  its  weapons  have  survived  prior  strikes. 

H.  Losses  because  of  unreliability  occur  before  the  bombers  reach  the  interceptors. 

These  assumptions  understood,  we  now  present  the  interceptor  defense  penetration 
model.  The  significant  parameter  is  PENPROB^y,  the  probability  that  an  individual 

independently  targetable  weapon  in  the  remaining  inventory  penetrates  the  interceptor 
defenses . 


Initially,  there  are  NNWP^y  bombers  of  type  IW  in  the  attacker’s  inventory.  At 
present,  there  are  NLEFT^y  bombers  that  have  not  been  sent  in  prior  strikes . The 

probability  that  these  NLEFT^  bombers  have  survived  prior  strikes  by  the  other  side 
is  SURV„„  . The  conditions  existing  on  the  three  strikes  are: 


IW 
NLEFT. 

SURV, 


!W  NNWPIW 


IW 
NLEFT 


= 1 


IW  " NNWPIW 


survt 


IW 
NLEFT 


<r  1 


IW  NNWPIW(i 


FSENTIW,  l> 


SURV 


IW 


s:  1 


KS  = 1 


KS  = 2 


KS  = 3 


(19) 


The  reprogrammable  reliability  (the  launch  reliability  or  the  probability  of  failing 
such  that  remaining  weapons  can  be  retargeted,  taking  the  failure  into  account)  for 
bomber  type  IW  on  stride  KS  is  RR^  ^ • The  non -reprogrammable  reliability  (in- 
flight reliability  or  the  probability  of  not  failing  after  it  is  too  late  to  retarget,  taking 
the  failure  into  account  for  bomber  type  IW  on  strike  KS  is  RN^y  ^g  . The  fraction  of 

the  remaining  bombers  of  type  IW  that  are  sent  on  strike  KS  is  FSENT^  ^g  . The. 
number  of  bomber  decoys  per  bomber  for  bomber  type  IW  is  BMRDECj^  . 

The  number  of  bombers  and  decoys  of  type  IW  successfully  launched  and  arriving 
at  the  defender’s  fighter  defense  system  on  strike  KS  is: 
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ARRIW  ' NLEFTIW*  SURVIW  ‘ fsentiw,ks  * rriw,ks 
. (1  + BMRDEC^) 


RNIW,KS 


(20) 


The  total  number  of  bombers  and  decoys  arriving  at  the  fighter  defense  system  is, 
then: 

B=£aRR  (21) 

IW  IW 

Now,  the  defense  includes  NIjp  fighters  of  type  IF  . It  is  assumed  that  a fraction  F 

of  the  bombers  is  to  attack  through  a given  corridor  and  that  the  same  fraction  F of  the 
fighters  is  pre -positioned  to  counter  the  F*B  bombers  in  the  corridor.  If  the  probability 
that  an  attacking  fighter  of  type  IF  will  kill  a bomber  of  type  IB^  is  PK^p  ^ , the 

probability  that  a specific  one  of  the  F • NIjp  fighters  encounters  and  kills  a specific 

bomber  out  of  the  F*  B in  the  corridor  is  PK.„  TR  / (F . B) 

» IW 


The  probability  that  none  of  the  F*  NI.p  fighters  kills  a given  bomber  of  this  type 

is,  then:  r i „ MI 

I F • NI 

IBj^ylF  = 1 ‘ 


PSURV 


TF 


(22) 


be  R 


Let  the  ratio  of  interceptors  of  type  IP  to  the  total  number  of  bombers  in  the  corridor 
IF  ‘ 


Rjp  = f-niif/(f-b) 


(23) 


The  probability  of  survival  is,  then: 
PSURV 


miwIF  = 


1 " (pkif,ib. 


IW 


/(F-B)) 


F-  B 


R 


IF 


(24) 


As  an  approximation  of  (24)  as  (F  • B)  becomes  large: 


R 


IF 


PSURV jg^jp  - [exp(-  PKif>ib^)] 

psur-Vib^if  = exP  C~  pi4f,  ib^  * niif/b] 


(25) 


When  all  interceptor  types  are  considered,  the  probability  that  a given  bomber  or  decoy 
penetrates  the  defenses  is  of  the  form: 

p 

PSURVm  = n exp  [ - A TF  JR  • (NIIp/B)  ’>IF’  IBIW  ] (26) 

KIW  IF  IF,1BIW  IF 
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Although  equation  (26)  represents  the  survival  probability  of  a bomber  or  decoy,  the 
important  issue  is  whether  a bomber  successfully  delivered  its  weapon  and,  if  so,  the 
reliability  of  those  weapons  after  release.  These  parameters  are  taken  into  account  in 
calculations  of  the  overall  penetration  probability  (PENPROBp^)  of  the  weapon  type. 


From  equation  (26),  the  probability  that  a bomber  or  decoy  of  type  IBp^  is  killed  is 


PKILLTn  = 1 - PSURVTn 
IBIW  BIW 


(27) 


The  fraction  of  bombers  of  weapon  type  IW  killed  before  weapon  delivery  is  FKBDj^, 
The  effective  probability  of  survival  to  weapon  launch  is,  then: 


PSEFF^  = PSURVffi[w  +(PKILlID[w  • (1  - FKBD^,)) 


(28) 


Finally,  tor  ASMs  the  probability  that  the  weapon  is  successfully  delivered  must 
take  the  ASM  reliability  (ARELp^)  into  account  as  a multiplication  factor.  The  overall 

penetration  probability  for  weapon  type  IW  is,  then: 


PENPROBM  = PSEFF^  • AREL^ 


’ ARELm  [l  -(fkbd^u  - psurv1B[w))] 


(29) 


The  number  of  bombers  and  decoys  of  weapon  type  IW  expected  to  penetrate  is: 


BDPEN^  = ARR^  PENPROB^ 


(30) 


However,  since  only  the  fraction  l/(BMRDECp^ 
expected  number  of  bombers  penetrating  is: 

PEN^  = ARR^  • PENPROBrw/(l  + BMRDEC^) 


+ 1)  are  actually  bombers,  the 

(31) 


Terminal  Defense  Penetration  and  Total  Price  to  Attack  a Target 

The  assumptions  regarding  surface -to -air -missile  (SAM)  defenses  are  as  follows: 


A.  Gravity  bombers  do  not  have  to  pay  a definite  price  to  attack  a SAM -defended 
target;  their  kill  probability  against  that  target  takes  into  account,  as  a multiplicative 
factor,  the  probability  that  the  target  will  shoot  down  the  bomber  with  SAMs. 


B.  ASM  bombers  with  short-range  attack  missiles  (SRAMs)  must  destroy  with  a 
probability  of  .95  the  requisite  number  of  suppressible  SAM  sites  and  exhaust  the 
requisite  number  of  missiles  from  non-suppressible  SAM  sites  before  attacking  the 
target. 


C.  Subsonic  cruise  armed  decoys  (SCADs)  must  exhaust  the  requisite  number  of 
suppressible  and  non-suppressible  SAM  sites  before  attacking  the  target. 
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D.  Each  SAM  may  kill  only  one  attacking  SRAM  or  SCAD. 

E.  Each  SRAM  may  attack  only  a single  suppressible  SAM  site. 

Based  on  the  above  assumptions,  gravity  bombers  have  no  entry  price  to  pay  at  a 
SAM-defended  target.  Therefore: 

NPjw  = 0 : gravity  bombers  (32) 

For  the  ASM  bomber  case,  the  price  for  suppressing  suppressible  sites  is  calculated 
first.  A single -shot  kill  probability  (SSK)  is  calculated  for  the  SAM  site  hardness,  the 
SRAM  yield,  CEP,  and  height  of  burst.  The  probability  of  kill  must  then  be  degraded  by 
the  probability  that  the  SRAMs  arrive.  The  resulting  overall  kill  probability  is,  then: 

SSKK  = SSK  • PENPROBRV  • RN^  Ks  (33) 

For  this  SSKK  value,  the  number  of  warheads  required  to  produce  at  least  a 95  per- 
cent kill  probability  (IWP^y)  is  obtained  from  the  relationship: 

.95  s 1 - (1  - SSKK)  m (34) 

Therefore: 

IWPj^y  = Iy  [log  (.05)/log  (1  - SSKK)]  (35) 

Now,  since  there  are  NSAMD^  suppressible  SAM  sites  at  target  IT,  the  total  number 
of  warheads  required  to  suppress  the  sites  is  IWP^,  • NSAMD. 

This  result  must  now  be  translated  into  independently  targetable  weapons  required 
per  target.  We  may  do  this  by  dividing  through  by  the  number  of  warheads  per  bomber 
of  weapon  type  IW  and  multiplying  by  the  number  of  independently  targetable  weapons  of 
weapon  type  IW  per  bomber,  and  taking  this  to  the  next  larger  integer  value,  to  obtain 
the  price  of  independently  targetable  weapons  for  suppressing  SAM  sites: 

NPSS^  = IyClWP^  * NSAMD  • TPRWPN^/WHDS^  ] (36) 

Next,  the  price  for  non-suppressible  SAM  sites  must  be  calculated.  This  is  very 
similar  to  the  missile  price  calculation,  since,  in  both  cases,  the  specified  number  of 
interceptors  must  be  exhausted. 

The  number  of  ASM  warheads  per  bomber  for  weapon  type  IW  is  WHDS^y.  The 
number  of  \SM  decoys  per  ASM  warhead  is  ASMDEC^  . The  total  number  of  warheads 
and  decoys  per  ASM  bomber  of  weapon  type  IW  is,  then,  WHDS^y(l  + ASMDEC^y) . 

Then,  since  the  NSAMX  non-suppressible  SAMs  are  capable  of  destroying  a single  decoy 
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or  warhead  each,  all  NSAMX  missiles  are  capable  of  destroying  BDES^  bombers  of 
weapon  type  IW: 

BDES^  = NSAMX/CWHDS^U  + ASMDEC^)]  (37) 

Taking  into  account  the  penetration  probability,  the  price  in  bombers  of  type  IW 
arriving  at  the  interceptor  defense  is: 

ARR*^  = NSAMX/ [PENPROB^  • WHDS^U  + ASMDEC^)]  (38) 

In  terms  of  bombers  successfully  launched,  thus  taking  inflight  reliability  into  con- 
sideration, the  price  is: 

ampm'arr^rniw,ks  <39> 

Then,  as  in  the  case  of  missiles,  if  this  result  is  multiplied  by  the  number  of  inde- 
pendently targetable  weapons  per  bomber,  the  result  is  expressed  in  terms  of  indepen- 
dently targetable  weapons  successfully  launched: 

NPSNIW  - \,[NSAMX  • TPRWPNlw/(RNlw  _ Kg  . PENPROBm  ^ 

•WHDSrw(l+ASMDECm)]) 


The  total  price  for  ASM  bombers,  then,  is  just  the  sum  of  the  prices  against  sup- 
pressible  and  non-suppressible  SAM  sites: 

NPW  = NPSS^  + NPSNpyy  (41) 

For  the  SCAD  case,  since  all  SAMs  must  be  exhausted  and  there  are  NSAMDM 
suppressible  SAMs  per  site,  the  total  number  of  SAMs  to  be  exhausted  is  (NSAMD 
• NSAMDM  + NSAMX).  Additionally,  for  the  SCAD  there  are  no  secondary  decoys,  and 
ASMDECj^  = 0.  The  overall  price  for  SCADs,  thus,  is: 

NP^y  = Iy  (NSAMD  • NSAMDM  + NSAMX)  . TPRWPN^ARN^  Kg 

• PENPROB^  • WHDS^)]  (42) 

The  number  of  independently  targetable  weapons  available,  NWP^  , is  just  the 

integer  number  of  bombers  that  will  be  successfully  launched  on  this  strike  (the  number 
left,  multiplied  by  the  fraction  sent  and  the  launch  reliability),  multiplied  by  the  number 
of  independently  targetable  weapons  per  bomber,  with  the  results  truncated  to  an  integer 
value: 

N^IW  = id[,d[nleftiw  • SURVIW  • fsentiw,ks  • rriw,ks!' TPRWPNIW 
L (43) 
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Now,  if  NP^  sNWP^y  , the  price  of  entry  exceeds  or  equals  the  available  resources, 
and  no  payoff  from  the  target  can  be  obtained  by  weapon  type  IW  alone.  If  NP^  < NWP^  , 
the  entry  price  allows  a payoff  from  weapon  type  IW. 

Sample -Bomber  Penetration  Problem 


A series  of  cases  will  now  be  considered  and  the  necessary  calculations  carried  out 
in  a tabular  manner,  with  one  case  explained  in  detail.  The  results  will  be  presented 
both  in  a table  and  graphically.  The  parameter  values  for  the  cases  are  as  follows: 


KS  > 1 
NLEFT 


IW  SURVIW 


rriw,ks  = 

.99 

RN  = 

*‘tw,ks 

.95 

BMRDEC^ 

= 1 

IBIW  = 

l 

WHDS^  = 

5 

TPRWPN^ 

= 5 

ASMDEC^ 

= 1 

NI1 

= 50 

Al,l 

= .5 

Cl,l 

= 1 

FKBD^ 

= .9 

= .9 

SSK 

= .5 

NS  A MX 

= 5 

NSAMD 

= 2 

NSAMDM 

= 5 

FSENT. 


IW.KS  = 10,20,30,40,50 


The  case  of  NLEFT^,  • SURV^  • FSENT^,  Kg  = 50  will  now  be  carried  out  in  detail. 
The  expected  number  of  bombers  and  decoys  arriving  at  the  interceptor  defenses  is: 


f i 
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ARRIW  = <NLEFTIW  ‘ SURVIW  * FSENTr\V,KS)  ’ rriw,ks  * rniw,ks 
• (1  + BMRDECrv) 

= (50)  • (.99)  • (.95)  • (2)  = 94.05 

Because  only  one  weapon  type  is  considered,  the  total  number  of  bombers  and  decoys 
arriving  at  the  interceptor  defenses  is: 

B = 94.05 

The  probability  that  a decoy  or  bomber  survives  the  interceptor  attack  is: 

PSURV1  = exp[-  .5  • (50/94. 05)1] 

= .7666 

The  penetration  probability,  then,  is: 

PENPROB^  = (.90)  • [l  - ((.9)  • (1  - .7666))] 

= .7109 

The  expected  number  of  bombers  penetrating  the  interceptor  defense,  then,  is: 
PEN^  = ((94.05).  (.  7109))/(2) 

= 33.43 

Now,  the  SAM  defense  must  be  handled.  First,  against  the  suppress ible  sites,  the 
overall  kill  probability  is: 

SSKK  = (.5)  • (.7109)  • (.95) 

= .3377 

The  number  of  ASMs  required  per  suppressible  site  is: 

WPiw  = Iy  [log  (.05)/log  (1  - .3377)] 

= ly  [7.27]=  8 

The  price  to  overcome  suppressible  SAM  sites,  then,  is: 

NPSSjw  = Iy  [(8).  (2)(5)/5  ] 

= Iy  [16]=  16 
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The  price  to  exhaust  the  non-suppressible  SAMs  is: 

NPSN^  = lu  j<5)  • (5)/ [(.95)  • (.7109)  • (5)  . (2)  ]} 

= I(J  [3.70]  = 4 

The  total  entry  price  to  the  target,  then,  is: 

NP^  = 16  + 4 = 20 

Now,  the  number  of  allocatable  weapons  was: 

NWPIw  = 'd  ['d  [<nleftiw • SURVIW  • fsentiw,ks>  • rriw,ks1  • ™W!,Niw 
• TPRWPNIWJ 

= ID[lD[(50)  . (.99)].  (5)J 
= 245 

Since  NP^  s NWP^,  a payoff  may  be  secured  from  the  target.  The  results  for  the 
other  cases  are  shown  in  table  2. 


TABLE  2 

SUMMARY  OF  RESULTS  FOR  SAMPLE  CASES 
IN  BOMBERS  UNLESS  OTHERWISE  NOTED 


NWP,W 
Independently 
targe table 
weapons 

ARR|yy 

PENPROB,w 

PEN  ,w 

SSKK 

IWP 

Independently 

targetable 

weapons 

NPSS,W 

Independently 

targetable 

weapons 

NPSN|W 

Independently 

targetable 

weapons 

NP,W 

Independently 

targetable 

weapons 

45 

18.81 

.2647 

.3044 

2.86 

.1446 

20 

40 

9 

49* 

95 

37.62 

.5145 

.5068 

9.53 

.2407 

11 

22 

6 

28 

145 

56.44 

.5421 

.6101 

17.22 

.2898 

9 

18 

5 

23 

195 

75.24 

.7173 

.6710 

25.24 

.3187 

8 

16 

4 

20 

245 

94.05 

.7666 

.7109 

33.43 

.3377 

8 

i 

16 

4 

20 

•No  payoff  from  this  caie,  iincaNP(w>  NWPjyy 
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SINGLE -SHOT  KILL  PROBABILITY  MODEL 


The  single -shot  kill  probability  model  has  several  distinct  parts.  The  first  involves 
determination  of  a 50  percent  single -shot  kill  radius  for  weapon  type  IT  when  used 
against  target  type  IW  . The  second  part  involves  the  introduction  of  aiming  errors  and, 
for  area  targets,  distribution  of  the  target,  to  determine  a probability  of  kill  for  weapon 
type  IW  against  target  type  IT  , on  the  premise  that  the  weapon  type  is  successfully 
delivered.  The  third  part  involves  strategy  considerations  and  the  probability  that  weapon 
type  IW  can  be  delivered  against  target  type  IT  . 

Symbols  used  in  the  equations  are  listed  below: 


A 

CD 

CEP 


F 

FM 

FM' 

FM” 

FMO 

F(t) 

GZ 

H 

HOB 

IMP 

K 

KN 

M 

MU 


PAT 

PKILL 


IW 


PK(R) 

PSURVIW,IT 


PT 


Parameter  for  generalized  kill  function 
Drag  coefficient 

Circular  error  probable  - radius  of  the  circle  centered  at  the  mean 
and  containing  50  percent  of  the  impact  points 

Adjustment  factor;  ratio  of  FM  for  any  yield  to  FMO  for  the  datum  yield 

Initial  value  of  the  force  applied  to  the  target 

Peak  force  per  unit  area 

Force  per  unit  area  of  blast  wave 

Peak  damage  pressure  for  a given  level  of  damage 

Force  resulting  from  dynamic  pressure  and  as  a function  of  time 

Ground  zero 

Lever  arm 

Height  of  burst 

Impulse  of  the  system 

Target  factor  used  in  reference  (c) 

Similar  to  the  K factor  in  reference  (c),  except  for  a factor  of  10 
Effective  mass  of  the  system 
Ductility  ratio 

Peak  overpressure  associated  with  a .5  probability  of  obtaining  a 
desired  degree  of  damage 

Atmospheric  pressure  ahead  of  blast  wave 

Probability  of  target  being  killed 

Probability  of  damage  as  a function  of  distance  R 
Probability  of  surviving  attack  by  the  target  terminal  defense 

Point  at  which  effective  mass  is  assumed  to  be  concentrated 
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P(t) 

Overpressure 

ptgtiw,it 

Probability  that  there  is  a surviving  target  at  the  site 

I i 

Q0.5 

Peak  dynamic  pressure  associated  with  a .5  probability  of  damage 

1 ‘ 

Q(t) 

Dynamic  pressure 

6- 

r 

R 

Distance  (see  figure  7) 

1 

R’ 

Distance  to  the  target  element  from  (0, 0) 

I 

R" 

Distance  to  the  impact  point  from  (0, 0) 

RE 

Structural  resistance 

§_ 

R(x) 

Resistance  of  target  to  displacement  surface  area  exposed  to  force 

s. 

S 

Surface  area  exposed  to  blast 

sskp1iw,it 

Single -shot  kill  probability  for  a single  warhead 

> 

i 

STRATjW  ^Strategy  factor 

p 

V 

1 

/ 

survx[T 

Probability  that  target  type  IT  has  not  already  been  destroyed 

1 

i 

r: 

T 

Natural  period  of  elastic  vibration 

l 

TC 

Period  of  a compound  pendulum 

* 

TD 

Effective  time  duration 

B 

w 

v 

TDO 

Effective  time  duration  for  a reference  yield  YO 

E 

TS 

Artificial  "period"  for  system 

T- 

V 

Velocity  of  the  system 

f 

A 

VF 

Maximum  and  final  velocity  of  system 

%- 

'v 

£ 

VN 

Vulnerability  number 

VO 

Initial  velocity  of  the  system 

■jr 

W 

Weight  of  the  system 

,U 

If 

WB 

Work  done  by  blast  wave  in  moving  the  structure 

* 

WR 

Weapon  radius 

Sr 

>- 

r 

ws 

Work  done  by  the  system  in  resisting  deflection 

c 

gr 

X 

Target  displacement  or  deflection 

1 

jT 

XE 

Yield  deflection 

|_ 

XM 

Maximum  deflection 

r_ 

Y 

Weapon  yield 

5 

YO 

Datum  yield 

a Angular  acceleration 

r)  Coefficient  of  friction  between  the  sliding  surfaces 

0 Kinetic  energy  of  the  system 

0’  Work  done  to  overturn  the  system 

(0  Angular  velocity 

(T  Measure  of  gradual  fall -off  in  probability  of  damage  with  distance 

or',  a"  Standard  deviations  of  the  distributions 

t Torque 

Determination  of  the  50  Percent  Kill  Radius* 

Blast  Phenomena  - The  blast  from  a nuclear  explosion,  which  constitutes  the  foremost 
cause  of  damage  to  structural  targets,  is  the  result  of  two  related  blast  phenomena:  over- 
pressure, P,  and  dynamic  pressure,  Q.  Overpressure  may  damage  targets  by  several 
different  modes:  by  crushing,  by  causing  severe  horizontal  deflections,  by  overturning, 
etc.  By  contrast,  dynamic  pressure  acts  in  one  direction  and  can  be  thought  of  as  wind 
pressure.  It  results  from  the  velocity  and  density  of  the  strong  transient  winds  that 
accompany  the  blast  wave.  Both  Pand  Q are  usually  expressed  as  pressure,  in  units  of 
pounds  per  square  inch  (PS1). 

When  the  blast  wave  arrives  at  a target,  both  P and  Q rise  abruptly  to  their  maxi- 
mum values  and  then  decay  as  functions  of  time,  in  the  manner  shown  by  the  solid  lines 
in  figure  1. 

It  is  more  convenient  to  work  with  pressures  which  vary  linearly  with  time  instead 
of  the  shapes  shown  by  the  solid  lines  in  figure  1.  The  linear  method  of  representing 
the  pressure -time  curves  uses  triangles  with  the  same  maximum  pressure  value  as  in 
the  ideal  case,  but  lasting  for  a new  duration,  termed  effective  duration,  TD  . This 
representation  is  indicated  by  the  dotted  lines  in  figure  1.  The  effective  duration  is 
usually  so  determined  as  to  provide  the  same  impulse  in  the  triangular  representation 
as  contained  in  the  ideal  case;  that  is,  the  areas  of  triangles  P-TD-0  and  Q-TD-0  are 
equal  to  the  areas  under  the  respective  ideal  curves.  Pressures  represented  in  this 
triangular  manner  are  referred  to  as  "initially  peaked  triangular  pressure  pulses." 

The  two  blast  phenomena,  overpressure  and  dynamic  pressure,  cause  damage  to 
targets  by  several  different  types  of  loadings.  Thus,  most  targets  may  be  classified 
into  two  types,  overpressure -sensitive  (or  P-type)  targets  and  dynamic -pressure -sensitive 
(or  Q -sensitive)  targets. 

Externally  Applied  Forces  - Although  many  targets  are  influenced  by  more  than  one 
form  of  blast  loading,  it  is  convenient  to  consider  targets  to  be  solely  either  overpressure 
or  dynamic -pressure -sensitive.  Also,  the  complex  decaying  pressures  that  result  in 
complex  forces  can  be  approximated  with  sufficient  accuracy  by  the  triangular  representa- 
tion discussed  previously  if  the  target  response  time  is  long  compared  with  the  duration  of 


* Except  as  noted,  this  derivation  is  based  on  reference  (a). 
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the  loading.  Since  the  loaded  area  of  many  targets  remains  constant  during  the  significant 
damage  phase,  these  pressures  are  directly  proportional  to  forces  on  the  target.  There- 
fore, these  forces  may  be  represented  triangularly,  as  were  the  pressures. 


FIG.  1 : OVERPRESSURE  AND  DYNAMIC  PRESSURE  AS  FUNCTIONS  OF  TIME 

Figure  2 is  a simplified  force-time  diagram,  where  F(t)  is  the  force  (pressure 
times  unit  area),  varying  linearly  with  time,  and  FM  is  the  initial  value  of  the  force 
applied  to  the  target.  This  peaked  force  could  result  from  either  of  the  two  pressures 
discussed  previously. 


FIG.  2:  SIMPLIFIED  FORCE  TIME  RELATIONSHIP 
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As  mentioned  above,  the  forces  on  the  target  during  the  significant  damage  phase  are 
equal  to  the  product  of  the  particular  pressure  times  the  area  existing  at  that  time.  For 
example,  the  existence  of  frangible  sheathing  may  be  an  important  determinant  of  the  area 
loaded.  In  the  special  case  of  a loading  resulting  from  dynamic  pressure,  the  shape  of  the 
target  is  important.  In  general,  the  force  resulting  from  the  dynamic  pressure  is  given 
by 

F(t)  = S • Q(t)  * CD  (44) 

where  F(t)  is  force  at  time  t,  S is  the  presented  area,  Q(t)  is  the  dynamic  pressure  at 
time  t,  and  CD  is  the  drag  coefficient.  This  coefficient,  a property  primarily  of  target 
shape,  ranges  in  value  from  0.35  for  cylinders  to  2 for  unshielded  structural  members 
(I-beams,  etc.).  It  also  depends  on  Reynolds  Number  for  most  objects. 

Resistance  of  Target  - Complex  targets  resist  displacement  in  the  same  manner 
as  a beam  resists  bending,  or  a desk  resists  sliding,  or  a box  resists  overturning.  The 
relationship  between  the  resistance  R(X)  that  many  targets  offer  to  displacement,  and 
target  displacement  or  deflection,  X may  be  expressed  by  a bi -linear  resistance  curve 
as  shown  in  figure  3.  In  figure  3,  the  point  (RE,  XE)  is  the  yield  point.  Resistance  func- 
tions are  shown  for:  (a)  a strain-hardening  situation  in  which,  for  some  reason,  the 
resistance  of  the  structure  increases  beyond  the  yielding  point  of  (RE,  XE),  (b)  a con- 
ventionalized elasto -plastic  relationship  with  constant  RE  after  the  yield  deflection  XE  is 
reached,  and  (c)  an  unstable  situation  in  which  the  resistance  decreases  after  the  yield 
point.  Although  resistance  curves  take  many  forms,  as  shown  in  figure  3,  they  may  be 
approximated  by  an  "equivalent"  elasto -plastic  resistance -deflection  relation.  In  such 
equivalent  cases,  shown  in  figure  4,  fictitious  values  of  RE  and  XE  are  used  so  that  the 
area  under  the  equivalent  curve  equals  the  area  under  the  original  elasto -plastic  curve; 
thus  the  work  required  to  deflect  the  structure  is  the  same.  To  assist  in  relating  external 
force  and  resistance,  the  structural  resistance,  R(X),  should,  for  convenience,  be  de- 
fined in  the  same  units  as  F(t),  i.e.,  pounds  or  psi. 


FIG.  3:  RELATIONSHIP  BETWEEN  STRUCTURAL 
RESISTANCE  AND  DISPLACEMENT 
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XM 


XE  XE' 

Displacement  (or  deflection),  X 

FIG.  4:  RELATIONSHIP  BETWEEN  EQUIVALENT 
RESISTANCE  AND  DISPLACEMENT 

Figure  4 introduces  a new  term,  XM,  called  the  "maximum  deflection, " which  is 
the  deflection  at  which  a target  is  damaged  to  the  extent  desired.  For  example,  in 
overturning  a target,  XM  will  be  the  point  at  which  the  target  is  balanced  with  its  center 
of  gravity  over  the  point  of  rotation.  In  another  case,  XM  might  represent  the  distance 
a bridge  slides  before  falling  off  its  supports:  or  XM  might  be  the  deflection  of  the  upper 
end  of  a building  column  at  which  it  will  fall  by  its  own  weight  or  the  building  will  be 
unusable  to  a specified  degree. 

In  order  to  describe  XM  conveniently,  the  term  ductility  ratio,  MU,  is  used.  It  is 
defined  as  follows: 

MU  = XM/XE  (45) 

The  ductility  ratio  is  a measure  of  the  inelastic  action  required  to  damage  a target. 

It  is  convenient  also  to  describe  structures  by  their  natural  periods  of  elastic  vibra- 
tion, T.  For  a single -degree -of -freedom  system,  the  period  is  represented  by: 

T = 2n(M«XE/RE)1/2  (46) 

where  M is  the  mass  of  the  structure,  and  RE/XE  is  the  spring  constant  of  the  structure. 
When  an  equivalent  ela3to -plastic  resistance  deflection  relationship  is  considered  as 
shown  in  figure  4,  the  resulting  period  is  then  termed  the  "effective  period." 

Equation  of  Motion  of  a Single -Degree -of -Freedom  System  - The  mass  and  resistance 
of  a structure  may  be  represented  by  the  single  degree -of -freedom  system  shown  in 
figure  5. 

Newton's  equation  of  motion  for  this  single -degree -of -freedom  system  may  be 
expressed  as: 

F(t)  - R(x)  = Md2X/dt2 


(47) 


where  F(t)  is  the  externally  applied  force  as  a function  of  time,  R(x)  is  the  resistance 
of  the  system  as  a function  of  displacement,  M is  the  effective  mass  of  the  system  and 
2 2 

d X/dt  is  the  second  derivative  of  displacement  with  respect  to  time  (acceleration). 


FIG.  5:  A SINGLE-DEGREEOF-FREEDOM  SYSTEM 


The  displacement  of  this  model  can  correspond  to  the  displacement  of  some  point  on 
a structure,  for  example,  point  PT  on  the  roof  of  the  one-story  structure  shown  in 
figure  6.  The  effective  mass  is  assumed  to  be  concentrated  at  this  point  PT.  The 
resistance  of  the  structure  is  related  to  the  sum  of  the  resistances  of  the  individual 
components  (the  columns)  of  the  structure.  The  external  force,  acting  on  the  concen- 
trated mass,  may  be  an  initially  peaked  triangular  force  of  the  type  described  in  figure  6. 


Applied  force 
F(t) 


Resisting  force 
R(X) 


FIG.  6:  STRUCTURE  BEFORE  AND  DURING  THE  BLAST 


The  solution  to  equation  (47)  for  structures  represented  by  an  elasto -plastic 
resistance  (figure  4)  was  solved  by  reference  (b);  the  derivation  is  given  below. 
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Derivation  of  the  Damage  Pressure  Level  Equation  - The  damage  pressure  duration 
for  most  targets  has  a value  somewhere  between  zero  (impulse  load)  and  infinity  (steady- 
state  load).  The  approach  is  to  add  together  the  determined  peak  damage  pressure 
equation  for  each  of  these  two  duration  limits.  This  results  in  a single  equation  for  a 
damage  pressure  level  with  a duration  somewhere  between  zero  and  infinity.  The  mathe- 
matical technique  of  adding  the  two  equations  results  in  a solution  that  closely  agrees  with 
the  exact  solution. 

The  basic  concept  of  the  derivation  is  the  Law  of  Conservation  of  Energy,  that  is, 
the  work  done  by  the  external  forces  from  the  blast  wave  must  equal  the  energy  absorbed 
by  the  structure  in  the  process  of  deformation.  This  means  that  the  work  done  by  the 
external  forces  is  equal  to:  (1)  the  work  done  in  deflecting  the  target,  or  (2)  the  work 
done  on  the  system  resisting  the  deflection,  or  (3)  the  kinetic  energy  of  the  system  for 
impulse  loading. 

Using  these  facts,  the  peak  damage  pressure  level  equation  for  impulse  loading  and 
then  for  the  steady -state  loading  will  be  developed. 

For  the  impulse  loads,  let  the  kinetic  energy  of  the  system,  0 , equal  the  work  done 
by  resisting  deflection,  WS  . The  kinetic  energy  of  a system  is: 

Q=M-V2/2  (48) 

where  Q is  the  kinetic  energy,  M is  the  mass  per  unit  area  perpendicular  to  the  blast 
wave,  and  V is  the  velocity  of  the  system.  To  solve  the  equation,  the  velocity,  V,  is 
first  found.  To  do  this,  the  impulse,  IMP,  is  first  considered.  The  impulse  is  equal 
to  the  area  under  the  pressure -time  curve  (figure  2). 

IMP  = FM'  • TD/2  (49) 

In  equation  (49),  FM’  is  the  peak  force  per  unit  area,  and  TD  is  the  effective  time 
duration.  Impulse  is  also  equal  to  the  change  in  momentum  of  a system  for  impulse 
loading. 


IMP  = M A V = M(VF  - VO)  (50) 

Also,  where  M is  the  mass  per  unit  area,  VO  is  the  initial  velocity  of  the  system  (VO  = 0 
when  at  rest)  and  VF  is  the  maximum  and  final  velocity. 

Thus: 


FM’  • TD/2  = M(VF  - VO) 

(51) 

A V = FM’  * TD/2M 

(52) 

Substituting  this  in  equation  (48),  then: 
0 = (M/2)  • (FM'  • TD/2M)2 


(53) 
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Referring  to  figure  4:  The  work  done  in  resisting  the  external  force  to  the  maxi- 
mum deflection,  XM  is  equal  to  the  area  under  the  elasto -plastic  resistance  curve: 

WS  = RE  • XM  - (RE  • XE)/2  (54) 

According  to  equation  (47): 

RE  • XM  - RE  • XE/2  = (M/2)  • (FM’  • TD/2M)2  (55) 

2 

Dividing  both  sides  by  RE  and  rearranging: 

FM'/RE  = (2/TD)  • j(M  • XE/RE)  • [(2XM/XE)  - l]j-1/2  (56) 

Equation  (56)  contains  the  terms:  M,  XE,  and  RE.  Conveniently,  these  same  terms 
are  included  in  the  fundamental  period  of  vibration,  an  expression  that  may  be  used  to 
describe  a structure. 

Substituting  the  ductility  ratio  MU  from  equation  (45)  and  the  period  T from 
equation  (46)  into  equation  (56)  yields: 

FM'/RE  = (T/(tt  • TD)  • (2MU  - 1)1/2  (57) 

This  completes  the  first  portion  of  the  damage  pressure  level  equation  for  impulse 
loading. 

A steady -state  pressure  (or  force  per  unit  area)  means  that  the  pressure  does  not 
decay  with  time  and  that  the  pressure  is  of  the  same  intensity  at  zero  deflection  as  at 
maximum  deflection,  SM  . Thus,  the  work  done  by  the  blast  wave  in  moving  the  struc- 
ture a distance  XM  by  a force  FM"  is  equal  to: 

WB  = FM"  • XM  (59) 

where  WB  is  the  work  done  by  the  blast  wave  in  moving  the  structure  and  FM"  is  the 
force  per  unit  area.  As  stated  before: 

WB  = W 

Thus:  FM"  • XM  = RE  • XM  - RE  • XE/2  (60) 

or:  FM"/RE  = 1 - XE/(2XM)  (61) 

Substituting  the  ductility  ratio  MU  = XM/XE  in  equation  (61)  yields: 

FM’VRE  = 1 - 1/(2MU)  (62) 

This  completes  the  second  portion  of  peak  damage  pressure  level  equation  for  a loading 
of  infinite  duration. 
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As  stated  before,  the  sum  of  the  two  limiting  expressions  is  a good  approximate 
solution  over  the  whole  range  of  durations.  Thus: 

FM/RE  = FM/RE  + FM"/RE 

1/2  (63) 

- (T/nTD)  • (2  MU  - l)i/z  - 1/(2MU)  + 1 

In  usage,  equation  (63)  has  been  found  to  be  in  error  by  less  than  20  percent  in  most 
cases  and  less  than  32  percent  in  all  cases  when  compared  with  numerical  solutions  of 
equation  (47). 

In  order  to  use  this  expression,  some  datum  yield  YO  is  selected  for  which  the  peak 
force  per  unit  area,  FMO,  and  the  duration,  TDO,  are  known.  Then  equation  (63)  can 
be  rewritten  for  the  datum  yield  as  follows: 

F MO/RE  =(T/nTDO)*  (2  MU  - 1)1/2  - 1/(2MU)  + 1 (64) 

The  ratio  of  FM  (for  any  yield)  to  FMO  (for  the  datum  yield)  is  termed  "Adjustment 
factor."  This  ratio  may  be  either  less  than  or  greater  than  unity.  Thus: 

R = FM/FMO  (65) 

The  following  expression  for  R can  be  derived  by  substitution  of  equations  (63)  and  (64) 
into  equation  (65): 


n _ FM/RE 

R " FMO/RE 

1/2  <66> 

_ (T/nTD)  ■ (2  MU  - l)1' 1/(2  MU)  + 1 

(T/ttTDO)  • (2  MU  - 1)1/2-  1/(2MU)  + 1 

By  rearranging  and  simplifying,  the  influence  of  weapon  yield  in  terms  of  duration  is  given 
by  the  following  expression: 


R = 1 


(T/ttTDO)  • (2  MU  - 1)J 


• (1  - TDO/TD) 


Since  the  expression  within  the  brackets  in  equation  (67)  is  a function  of  the  target 
type's  physical  characteristics  and  the  effective  duration  for  the  datum  yield,  equation 
(67)  may  be  rewritten  as: 

R = 1 - KN  • (1  - TDO/TD)  (68) 

The  term,  KN,  is  similar  to  the  K factor  in  AFM  200-8,  except  for  a factor  of  10. 
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Thus: 


K = 10  KN  (69) 

Then,  equation  (68)  can  also  be  written  as: 

R = 1 - .IK  • (1  - TDO/TD)  (70) 


and: 


^ _ 10 
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(T/(n  • TDO))  • (2  MU -1) 1/2 

All  of  the  factors  included  in  K are  constant  after  the  datum  yield  is  selected. 

Since  much  of  the  original  target  damage  data  was  obtained  from  the  Japanese  explosions 
and  from  nominal  yield  experiments,  a datum  yield  of  20  kilotons  was  selected  for 
AFM  200-8. 

Hence,  KN  or  K is  a constant  for  a given  damage  level  to  a given  target.  This 
constant  K then  can  be  used  to  describe  the  way  in  which  target  damage  pressure  (FM) 
changes  as  a function  of  weapon  yield  with  respect  to  the  damage  pressure  level  at 
20  kilotons . 

The  K factor  above  is  for  a target  with  an  elasto -plastic  or  an  equivalent  elasto- 
plastic  resistance.  Different  K factor  relations  are  required  for  sliding  and  overturning 
targets  because  they  have  other  than  elasto-plastic  resistance.  These  relations  are 
derived  in  the  next  two  subsections. 

Derivation  of  K for  Overturning  Targets  - The  overturning  K factor  equation  is 
different  from  the  one  just  described  primarily  because  the  overturning  target  does  not 
follow  the  conventional  elasto-plastic  resistance  pattern.  This  subsection  will  develop 
the  overturning  K factor  equation  which  will  be  used  in  the  same  way  as  the  previously 
developed  K factor. 

Here  again,  the  approach  used  in  reference  (b)  is  to  determine  the  peak  damage 
pressure  equation  for  an  impulse  load  and  then  for  a steady -state  load.  These  two 
equations  are  added,  to  provide  a single  equation  for  a peak  damage  pressure  for  a 
duration  between  zero  and  infinity.  This  single  equation  will  be  used  to  determine  an 
adjustment  factor  relationship  after  which  the  overturning  target  K factor  may  be 
isolated. 

Conservation  of  energy  shows  that  the  work  done  by  an  external  force  from  the 
blast  wave  is  equal  to:  (1)  the  work  done  in  overturning  the  target  and  (2)  the  kinetic 
energy  of  the  system  for  impulse  loading. 

To  obtain  the  damage  pressure  level  equation  for  impulse  loading,  the  kinetic 
energy  equation  will  first  be  developed.  Then  the  kinetic  energy  will  be  equated  to  the 
work  done  in  overturning  the  target. 
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Consider  the  overturning  target  described  in  figure  7. 


The  following  assumptions  are  made:  (1)  the  target  overturns  about  point  0;  (2)  the 
target  does  not  slide,  and  (3)  the  duration  of  the  shock  is  short  compared  with  the  time 
required  for  the  target  to  overturn. 

From  energy  considerations,  the  following  kinetic  energy  equation  is  given  for  an 
impulse  load  for  an  initially  peaked  triangular  pulse  that  is  very  short  (impulsive)  in 
duration: 


KE  = 0=Io)2/2  (71) 

where  I is  the  polar  moment  of  inertia  and  o>is  the  angular  velocity  of  the  target.  For 
convenience: 

Q = (1  co)2/2I  (72) 

Another  relationship  from  physics  is: 

t = aT  (73) 

where  a is  the  angular  acceleration  of  the  target,  I the  polar  moment  of  inertia  of  the 
rotating  system,  and  t the  torque  (or  the  moment  equal  to  the  force  times  the  lever 
arm).  In  the  overturning  problem,  the  force,  or  force/unit  area,  i.e.  pressurt;,  is  a 
function  of  time  as  shown  in  figure  1 or  2.  Thus,  the  torque  in  equation  (73)  may  be 
expressed  in  terms  of  a moment: 

I • a = F(t)  * H (74) 

Since  the  integration  of  angular  acceleration  results  in  angular  velocity,  then  by  integra 
ting  both  sides,  equation  (74)  becomes: 

1 — / ,p»  • H • dt  (75) 
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The  pressure  or  force  at  any  time,  t,  between  the  limits  of  t = 0 and  t = TD,  is 
represented  in  figure  8 and  by  equation  (77). 


FIG.  8:  RELATIONSHIP  BETWEEN  FORCE  (PRESSURE)  AND  TIME 


FM/F(t)  = TD/(TD  - t) 


F(t)  = FM  * (TD  - t)/TD 


Substituting  the  force  (or  pressure)  at  any  time,  t,  between  the  limits  of  t = 0 and 
t = TD  into  equation  (75)  yields  the  result: 


-TD 

I • cO  = / FM  • H • [(TD  - t)  - /TD]  • dt 

o 

= FM  • H • t - FM  • H • t2/(2  • TD)  JD 


Consider  the  static  moments  about  point  0 in  figure  7.  The  resultant  static  force  acts 
at  a distance  H above  che  surface.  This  force  must  be  large  enough  to  overcome  the 
static  resistance  in  order  to  overturn  the  target. 


Thus: 


RE  . H = W • XC 


where  RE  is  considered  the  static  force  required  for  overturning,  (or  the  internal 
resistance),  H is  the  point  of  application,  and  W is  weight  of  the  target.  Substituting 
equation  (79)  into  (78),  and  modifying  FM  to  FM'  to  indicate  a very  short  pulse,  the 
equations  become: 
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(80) 


I • CO  = (FM’/RE)  • RE  • H • TD/2 

= (FM’/RE)  * (W  • XC/H)  • (H  • TD/2) 

Substituting  equation  (80)  into  (72),  the  solution  is: 

0 = (l/(2  • I))  • (FM’/RE)2  . (W  • XC  • TD/2)2  (81) 

With  the  kinetic  energy  solved,  the  work  done  in  overturning  must  be  determined.  The 
work  done  on  the  system  to  overturn  it  is  equal  to  lifting  its  weight  from  the  position  of 
rest  to  a position  such  that  the  center  of  gravity  is  directly  above  the  center  of  rotation 
as  shown  in  figure  9. 
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FIG.  9:  POSITIONAL  RELATIONSHIPS  FOR  AN  OVERTURNING  TARGET 


[D2~(XC)2]’/j 


0 


'=w{d-[D2-(XC)2]1/2} 

= W • D • |l  - [1  - (XC/D)2  ]1/2  j 


(82) 


Expanding  the  square  root  in  binomial  series  and  dropping  the  higher  order  terms: 

O’  = W • D • |l  - 1 + (XC/D)2/2  | 

= W • D * (XC/D)2/2  (83) 

= W • (XC)2/(2  • D) 


As  stated  previously,  fi  from  equation  (81)  and  Q'  from  equation  (83)  may  be 
equated: 

(1/21)  • (FM’/RE)2  • (W  • XC  • TD/2)2  = W • (XC)2/(2  • D)  (84) 

Therefore: 

FM'/RE  = 2 • (1/(D  • W))1/2/TD 


(85) 
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When  the  duration  of  the  force  pulse  is  very  long,  the  peak  force  must  at  least 
equal  the  initial  resistance  (RE)  or  tlte  system  will  not  overturn.  Thus  FM"  = RE: 

FM"/RE  = 1 (86) 

where  FM"  indicates  a force  of  great  duration. 

The  damage  pressure  for  an  initially  peaked  triangular  force  pulse  with  a duration 
between  that  of  an  impulse  (TD  = 0)  and  infinity  (TD  = 0)  can  be  approximated  by  the  sum 
of  equations  (85)  and  (86). 

FM/RE  = (FM'/RE)  + (FM”/RE) 

1/2  <87> 

= 1 + 2*(i/(D  • W))1//TD 

To  arrange  the  solution  in  terms  of  a K factor,  consider  the  period  of  an  over- 
turning object.  This  period  corresponds  to  that  of  a compound  pendulum  supported  at 
the  point  about  which  it  pivots  when  it  overturns. 

A compound  pendulum  (or  physical  pendulum)  is  any  body  that  vibrates  in  a 
manner  of  a pendulum  but  whose  mass  is  distributed  throughout  the  body,  and  not  a 
concentrated  mass  at  the  end  of  a cord  of  negligible  weight,  as  is  the  case  with  the 
simple  pendulum.  Thus,  the  period  of  a compound  pendulum  is: 


since: 


Then: 


TC  = 2tt  • (I/(M  • 


M = W/g  . 


D-g)) 


TC/2tt=  (I/(W  • D))1/2  (90) 

Conveniently,  the  terms  necessary  to  describe  an  overturning  target  in  the  period  ex- 
pression are  also  in  the  force  calculations  for  overturning.  Therefore  by  substituting 
equation  (90)  into  (87),  the  result  is: 

FM/RE  = 1 + TC/(tt  • TD)  (91) 

The  adjustment  factor  (which  leads  to  tiie  K factor)  for  overturning  targets  may  be 
solved  by  substituting  equation  (91)  into  the  previously  described  adjustment  factor  ex- 
pression, equation  (65). 

F = FM/FMO 
= (FM/RE)/(FMO/RE) 

_ 1 + TC/(tt  . TD)  (92) 
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By  rearranging  and  simplifying: 

F - 1 - [ i-r^TWTC)] - (1  - TD0/TD)  <*» 

Since  the  expression  within  the  brackets  in  equation  (93)  is  a function  of  the  target 
type's  physical  characteristics  and  the  effective  duration  for  the  datum  yield,  equation 
(93)  can  be  rewritten  in  the  form  of  equation  (68), 

F = 1 - KN  ♦ (1  - TDO/TD)  (94) 

By  again  making  K = 10  • KN,  equation  (94)  becomes: 

F = 1 - .IK  • (1  - TDO/TD)  (95) 

K = 10/ [ 1 + (rr  • TDO/TC)  ] 

Derivative  of  K for  Sliding  Targets  - Sliding  targets  have  a different  type  K factor 
expression  from  the  one  described  previously,  because  the  sliding  system  has  a rigid- 
plastic  resistance  (figure  10).  To  develop  the  sliding  target  K factor,  the  approach  used 
will  be  the  same  as  in  for  overturning  targets,  which  is  to  determine  the  peak  damage 
pressure  equation  for  an  impulse  load  and  then  for  a steady -state  load.  These  two 
equations  are  added  together  to  provide  a single  equation  for  a peak  damage  pressure 
for  a duration  between  zero  and  infinity.  This  single  equation  will  be  used  to  determine 
an  adjustment  factor  relationship  after  which  the  sliding  target  K factor  may  be  isolated. 

When  considering  sliding  as  the  mode  of  damage,  the  target  will  not  begin  to  slide 
(deflection  remains  zero)  until  the  applied  external  force  becomes  greater  than  the 
sliding  resistance  which  is  presumed  to  remain  constant.  An  example  of  a sliding 
target  would  be  a bridge,  which  may  alide  a maximum  distance  of  XM  on  its  abutment. 


Displacement  {or  deflection),  X 


FIG.  10:  RELATIONSHIP  BETWEEN  RESISTANCE  AND 
DISPLACEMENT  FOR  A SLIDING  TARGET 
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The  energy  absorbed  by  the  system  is  equal  the  area  under  the  resistance  curve,  or: 

0 = RE  • XM  (96) 

If  the  initially  peaked  triangular  force  pulse  is  of  very  short  duration  (impulsive) 
then  the  energy  absorbed,  0,  in  deformation  must  be  equal  to  the  impulsively  applied 
kinetic  energy,  O'  , or: 

0=0'  (97) 

Since  I = M • V for  impulsive  loadings: 

O'  = M • V2/2 

= l2/(2  • M)  (98) 

Also,  impulse  is  equal  to  the  area  under  the  initially  peaked  triangular  pressure- 
time curve,  or  considering  the  pressure  as  a force  per  unit  area,  then: 

1 = FM’  • TD/2  (99) 

Substituting  equation  (99)  into  equations  (97)  and  (98)  yields: 

0=0'  = (FM*  • TD/2)2/2M  (100) 


RE  • XM  • 2M  = (FM'  . TD/2)^ 
and  therefore: 

FM'/RE  = 2 • (2  • M • XM/RE)1/2/TD  (101) 

where  FM'  indicates  a very  short  duration  loading. 

When  the  duration  of  the  loading  is  very  long  (TD  -»  ®),  then  the  peak  pressure  or 
peak  force  (FM)  must  be  equal  to  or  greater  than  the  sliding  resistance;  otherwise,  the 
structure  will  not  move. 

Thus: 


FM’VRE  = 1 (102) 

According  to  reference  (b),  the  damage  pressure  level  for  an  initially  peaked  tri- 
angular force  pulse  with  a duration  between  that  of  an  impulse  and  infinity  can  be  approxi- 
mated by  the  sum  of  equations  (101)  and  (102),  or: 

FM/RE  = (FM'/RE)  + (FM"/RE) 

= 2 • (2  • M • XM/RE)1/2/TD  + 1 (103) 
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The  incorporation  of  this  expression  in  terms  of  the  previously  developed  adjustment 
factor  will  provide  a K factor  for  sliding  targets.  The  adjustment  factor  expression, 
equation  (65)  and  equation  (103)  may  be  written  as  follows: 


F = FM/FMO 
= (FM/RE)/(FMO/RE) 

_ 2 . (2  . M • XM/RE)1/2/TD  + 1 
2 . (2  • M • XM/RE)1/2/TDO+  1 


(104) 


As  shown  in  figure  10,  the  sliding  resistance  is  constant,  or  XE  can  be  taken  to  be 
zero.  Therefore  from  equations  for  the  ductility  ratio  (equation  (45))  and  natural  period 
of  vibration  (equation  (46))  T = 0 and  MU  = » . However,  to  incorporate  the  sliding 
target  into  the  K factor  system,  the  following  conveniently  artificial  "period"  may  be 
used: 


TS  = 2 • tt  . [2 . (W/g)  . (XM/RE)]1/2 
or 

TS/(2tt)  = [2*  (W/g)  • (XM/RE)]1/2 
Substituting  equation  (106)  into  equation  (104): 


t TD+  1 

7TD0  + 


1 


(105) 

(106) 


(107) 


Rearranging  and  simplifying  yields: 


F 


1 

1 + tt  • rDo/T$ 

J 


(1  - TDO/TD) 


(108) 


As  in  the  previous  cases,  equation  (108)  can  be  separated  into  two  parts,  one  of  which 
is  only  a function  of  the  target  type's  physical  characteristics  and  the  effective  duration 
for  the  datum  yield  and  can  be  rewritten  as: 


F = 1 - KN  • (1  - TDO/TD)  (109) 

For  K = 10KN  equation  (109)  becomes: 

F = 1 - .IK  • (1  " TDO/TD)  (110) 

K = 10/ [l  + (tt  .TDO/TS)] 


Considering  equation  (105)  once  again,  W is  the  weight  of  the  target,  XM  is  the  displace- 
ment of  the  target  and  RE  is  the  sliding  resistance  offered  by  the  target.  Assuming  the 
only  sliding  resistance  is  that  offered  by  friction,  then  RE  = r|W  where  q is  the  co- 
efficient of  friction  between  the  sliding  surfaces.  Using  this  assumption,  the  TS  term 
in  equation  (110)  becomes: 
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TS  = 2 • tt  • [(2  • W/g)  • (XM/(n  • W))]1/2 
= ' • tt  . (2  . XM/(g  • n))1/2 


(111) 


Effective  Time  Durations  and  VN  Number  Adjustments  - As  is  well  known,  a target's 
susceptibility  to  nuclear  weapon  blast  effects  can  be  described  in  terms  of  a vulnerability 
number,  consisting  of  three  parts:  a number,  a target  type  (i.e.,  P or  Q)  and  a K factor. 
The  latter  terms  have  been  described  previously  but  a cursory  review  of  the  vulnerability 
number  (VN)  scheme  is  in  order. 

The  numerical  value  scale  of  the  VN  is  an  arbitrary  classification  scheme  describ- 
ing a target's  vulnerability.  (The  basic  VN  is  referenced  to  20  kilotons  while  the 
"adjusted  VN"  is  adjusted  to  some  other  particular  yield.)  Actually,  the  numerical 
values  are  linear  functions  of  the  logarithm  of  either  the  peak  overpressure  or  the  peak 
dynamic  pressure  that  is  required  to  achieve  a given  probability  of  a defined  level  of 
damage  to  a randomly  orientated  target.  In  the  special  case  of  a .5  probability  of  damage, 
the  numerical  value  of  the  "P"  type  VN  is  defined  as: 


’P"VN  = 


log  (PQ>5)  - log  (1.1216) 
log  (1.2) 


(112) 


where  Pq  _ is  the  peak  overpressure  associated  with  a .5  probability  of  obtaining  a 

desired  degree  of  damage  such  as:  "severe,"  "light,"  "depot  repair,"  "30-day 
recovery  time,"  etc.  In  a similar  way,  the  "Q"  type  VN  is  defined  as: 


log  (Qn  s)  - log  (0.02893) 
"Q"VN  = 


(113) 


where  Qq  _ is  the  peak  dynamic  pressure  associated  with  a .5  probability  of  damage, 
These  two  relations  are  illustrated  in  figure  11. 


As  a point  of  interest,  it  may  be  noted  that  for  a given  numerical  value  of  either 

type  VN  (e.g.,  VN-12P  and  VN-12Q)  the  corresponding  overpressure  (Pn  . = 10  psi) 

u • o 

and  dynamic  pressure  (Qq  - = 2.3  psi)  are  related  by  an  approximation  given  by 
equation  (114) 


Q = 5 • P /(2  • (7  PAT  + P)) 
= .023  * P2 


(1 14) 


In  equation  (114)  PAT  is  the  atmospheric  pressure  ahead  of  the  blast  wave.  This 
means  that  under  nearly  ideal  surface  burst  conditions  and  with  no  K factor  considera- 
tions, in  the  example  cited,  VN-12P  and  VN-12Q  targets  will  be  damaged  with  a proba- 
bility of  .5  about  the  same  distance  from  ground  zero  (no  weapon  delivery  errors 
considered). 
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FIG.  11 : RELATIONSHIP  BETWEEN  50  PERCENT  PROBABILITY  OF 
PRESSURE  DAMAGE  AND  VULNERABILITY  NUMBER 


Occasionally,  the  .5  probability  of  damage -distance  is  used  erroneously  to  indicate 
a "go -no  go"  or  "cookie -cutter"  probability  of  damage -distance  situation.  Actually,  a 
target's  destruction  is  a "go  or  no-go"  proposition  only  when  considering  a given  damage 
criterion,  such  as  75  percent  building  floor  space  unusable  or  a vehicle  is  overturned. 
This  means  that  a target  is  either  damaged  or  not  damaged  to  the  desired  extent.  How- 
ever, achievement  of  this  criterion  must  be  considered  on  a probabilistic  basis  because 
of  variabilities  in:  weapon  output  and  blast  phenomena,  intrinsic  target  structural 
characteristics,  mutual  target  sliielding,  terrain,  target  orientation,  etc.  When  a 
specific  target  is  considered,  these  factors  lead  to  a variation  in  probability  of  damage 
versus  distance  from  ground  zero  (GZ),  as  shown  in  figure  12, 

It  may  be  noted  that  the  distance,  cr,  between  a probability  of  .69  and  .31  is  a 
measure  of  this  gradual  fall -off  in  probability  of  damage  with  distance.  Although  a .5 
probability  of  damage  distance  can  be  used  as  a "probable  radius,"  the  VN  system 
employs  another  distance,  which  is  the  weapon  radius  or  WR  . 


WR=  [J"  2 PK(R)RdR]1/2 


(115) 


where  PK(R)  is  probability  of  damage  as  a function  of  the  distance  R . The  WR  may  be 
thought  of  as  defining  the  radius  of  a circle  centered  on  the  GZ  . Assuming  an  infinite 
array  of  identical  targets,  the  WR  circle  will  enclose  just  as  many  undamaged  targets 
as  will  be  damaged  outside  the  circle.  Special  "cookie  cutter"  cases  of  a = 0 result 
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in  WR  equal  to  the  ’'probable  radius.”  For  real  conditions,  such  as  existed  at  Hiroshima 
and  Nagasaki,  <r/WR  may  very  well  be  0.20  or  0.30  ("Sigma-20"  or  "Sigma -30").  For 
Sigma-20,  the  solution  of  equation  (115)  leads  to  a WR  equal  to  a distance  R from  the 
GZ  corresponding  to  a probability  of  damage,  PK(R),  of  .473  . Thus,  for  many  real- 
life  situations,  WR  is  slightly  greater  than  the  .5  probability  radius. 


Distance  from  ground  zero 

FIG.  12:  EFFECT  OF  DISTANCE  FROM  GROUND  ZERO 
ON  PROBABILITY  OF  KILL 


Obviously,  the  described  probability  of  damage -distance  relationship  may  be  ex- 
pressed in  a probability  of  damage -pressure  relationship  shown  in  figure  13  for  either 
P or  Q,  by  converting  distance  into  pressure  by  empirical  overpressure  or  dynamic 
pressure -distance  relations. 
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A statistical  analysis  of  the  probability  of  damage -pressure  data  obtained  from 
Hiroshima  and  Nagasaki  was  performed  for  the  VN  system.  These  results  were  com- 
bined with  the  previously  described  VN-pressure  definitions  which  resulted  in  probability 
of  damage -VN-pressure  relationships. 

As  a matter  of  interest,  it  may  be  noted  that  in  reference  (c),  these  two  probability- 
pressure-VN  curves  were  "tied  together,"  that  is,  the  dynamic  pressure  and  overpressure 
were  related  at  all  points  by  the  approximate  form  of  the  Rankine -Hugoniot  equation 
2 

(Q  = .023P  ).  Originally,  only  the  overpressure  relation  was  based  upon  empirical 
evidence  and  the  dynamic  pressure  relation  merely  calculated  from  it.  A recent  restudy 
of  Hiroshima -Nagasaki  data  combined  with  DASA  field  test  data  resulted  in  confirmation 
of  the  probability -overpressure  curves  but  provided  new  probability -dynamic  pressure 
curves.  These  curves  are  interrelated  through  the  approximate  Rankine -Hugoniot  relation 
at  the  P and  Q _ points  only. 

Therefore,  a given  VN  represents  a line  on  a probability-pressure  graph  or  a curve 
on  a probability -distance  graph,  rather  than  a specific  point,  as  may  seem  to  have  been 
implied  in  equations  (112)  and  (113).  Thus,  while  reference  is  often  made  to  the  P or 

Q 5 corresponding  to  a VN,  pressures  associated  with  other  probabilities  of  damage  may 

be  read  from  probability  of  damage -VN-pressure  curves. 

Of  more  ultimate  value  is  the  use  of  these  curves  to  calculate  WRs  for  a wide  range 
of  weapon  height  of  burst-VN  combinations.  This  is  done  by  merely  retracing  the  steps 
described  above.  Specifically,  the  relationships  are  determined  at  one  height  of  burst 
(HOB)  and  are  assumed  to  remain  the  same  regardless  of  HOB.  Then,  P and  Q versus 
distance  data  for  one  kilotron,  as  determined  by  field  measurements  for  various  HOBs 
are  used  to  get  probability -versus -distance  curves  analogous  to  those  in  figure  12.  The 
integral  of  the  WR  equation  (114)  is  then  found.  In  practice,  this  has  been  accomplished 
by  electronic  computers  and  the  results  are  presented  in  graphs  which  are  shown  in  the 
forms  in  figure  14. 

These  one-kiloton  curves  can  be  used  for  any  nuclear  weapons  yield  by  the  proper 
scaling.  WR  scales  with  the  cube  root  of  the  weapon  yield,  so  that: 

WR  = (WR1KT)  (Yi/3)  (116) 

K factor  equations  using  the  approximate  effective  time  duration  for  both  dynamic 
pressure  and  overpressure  of  TD  = O.lYl/3  were  formerly  used.  The  use  of  this  simpli- 
fied relation  for  TD  was  entirely  satisfactory  for  targets  vulnerable  to  pressures  of  the 
order  of  10  psi  overpressure  and  for  the  pressure  data  available  at  that  time.  The  avail- 
ability of  more  recent  pressure  data  plus  the  increased  interest  in  hard  targets  as  well 
as  fractional  KT  nuclear  weapons  necessitated  the  improvement  of  the  approximate 
effective  duration  equation.  In  this  subsection  the  necessary  modification  of  the  effective 
duration  equations  for  both  overpressure  and  dynamic  pressure  will  be  explained. 
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Since  the  effective  duration  of  the  overpressure  and  the  dynamic  pressure  differ, 
separate  adjustment  factor  relationships  will  be  introduced  for  each.  Using  the  modified 
effective  durations,  revised  K factor  equations  will  be  provided. 


Weapon  radius 

FIG.  14:  EFFECT  Of  HEIGHT  OF  BURST  ON  WEAPON  RADIUS 

FOR  aone-kiloton  yield 
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Overpressure  Effective  Duration  and  Adjustment  Factor  - The  meaning  of  effective 
duration,  TD,  was  explained  in  the  subsection  on  determining  the  50  percent  kill  radius 
and  its  relationship  to  the  positive  phase  duration  shown  in  figure  1. 

Empirical  data  (reference  (d»  indicate  the  following  equation,  for  any  height  of  burst 
to  be  a better  approximation  for  a burst  of  yield  of  Y kilotons. 

TD  = .45  • p-1/2yl/3  (117) 

As  a result  of  equation  (1 17),  the  ratio  of  effective  duration  for  yield  Y to  the  effective 
duration  TDO  for  a reference  yield  YO  is  given  by: 

TD/TDO  = (PO/P)i/2  • (Y/YO)l/3  (118) 

The  adjustment  factor  for  damage  by  elasto -plastic  deformation,  overturning  and 
sliding  was  shown  to  be: 

T = 1 - . 1K(1  - TDO/TD) 

By  substituting  equation  (118)  into  equation  (119): 

F = 1 - , IK  [l  (P/PO)1/2  • (YO/Y)1/3] 

Now,  by  the  definition  of  the  adjustment  factor  R: 

R = FM/FMO 
= S • P/(S  • PO) 

Therefore,  equation  (120)  can  also  be  written: 

F = 1 - .IK  [1  - F1/2  • (YO/Y)1/3] 

1/2 

By  letting  X = F , equation  (122)  may  be  written  as: 

X2  - . IK  • (YO/Y)1/3  • X + . IK  - 1 = 0 (123) 

Because  the  values  of  K are  0 iK  s9,  equation  (123)  can  be  shown  to  have  two  roots, 
one  positive  and  one  negative.  Of  these  two,  only  the  positive  root  fits  the  physical  condi- 
tions . Therefore: 


(119) 

(120) 

(121) 

(122) 


, IK  • (YO/Y)1/3+  ([.IK  - _(YO/Y)1/3]2  + 4 (1  - , IK)) 1/2 


and  as  a result: 


(124) 


F = 1/4  j.  IK-  (YO/Y)1/3+  {[.IK  • (YO/Y)1/3]2  + 4(l  - ,1K)}1/2|2  (125) 
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duration3 TDQUfor  Preference  ^eld  ^O^i^given  by^^6  dUration  for  yield  Y t0  the  effective 


(128) 


TD/TDO  = (QO/Q)  1/3(Y/YO)1/3  (1 

The  equation  for  the  dynamic  pressure  adjustment  factor  becomes: 

F = 1 - . 1K(1  - TDO/TD) 

= 1 - .lK[l  - (Q/QO)1//3  . <YO/Y)1//3]  (128) 

Now,  because  of  the  definition  of  F: 

F = FM/FMO 

= S • Q/(S  • QO)  (129) 

Therefore,  by  using  the  ratio  of  Q/QO  from  equation  (128)  in  equation  (129): 

F = 1 - .1K[1  - F1/3  . (Y/YO)1/3  ] (130) 

By  letting  Z = F ^ , equation  (130)  may  be  written  as: 

Z3-.1K-  (YO/Y)1/3-  Z + . IK  - l = 0 (131) 

First  note  that  equation  (131)  has  only  a single  change  of  sign  since: 

- .1K(Y0/Y)1/3  <0 

.IK  - 1 <0  (132) 

° ne  used  to  find  (P>5)y  or  (Q>5)y  respectively  for  a weapon  of  yield  Y.  First,  an 

adjusted  VN  number  is  defined  for  the  P and  Q cases  respectively  as: 


tidbit » 


=E=.V  S3? ' 


* fspsa^5 


ADfP"VN  = 
ADJMQ*,VN  = 


log<P  -)Y  - log  (1.1216) 

T5TTJJ 

los  (Q.5)y  " log  (°*02893) 

log  (1.44) 


Rearranging  equations  (133)  yields: 


<P.5>Y  = (1.1216).  (i.2)ADr™ 

(Q  5)y  = (0.02893)  . (1.44)ADJ,,Q”VN 


(133) 


(134) 


Now,  by  the  definition  of  F : 


F = 


<P.5>Y 
^ . 5*  20  KT 


(1.1216)  . (1.2)ADr,P,,VN 
(1.1216)  .(1.2)"P‘VN 


(135) 


F = 


(Q.5^  20  KT 


(0.02893)  . (1.44)AD],,Q"VN 
(0.02893)  . (1.44)”^I,VN 


Solving  for  the  adjusted  VN  number: 


ADJ”P’VN  = "P"VN  + log  (Fj/Iog  (1.2) 
ADJ"Q"VN  = ”Q"VN  + log  (F)/log  (1.44) 


(136) 


Now,  tables  relating  weapon  radii,  WR  , to  VN  number  for  various  heights  of  burst 
based  on  a one-kiloton  yield  are  readily  available.*  The  radii  obtained  from  these  tables, 
however,  must  be  scaled  using  the  cube  root  scaling  laws  for  yield.  Thus,  the  weapon 
radius  for  yield  related  to  a one-kiloton  weapon  radius  is: 


WRy/WR1  = Y1/3 

Combining  equations  (136)  and  (137)  as  advocated  by  equation  (116): 
WRy("P"  CASE)  = Y1/3  • WR1  (BASED  ON  "P"  CASE  ADJ  VN) 

WRy("Q"  CASE)  = Y1/3  • WR1  (BASED  ON  "Q"  CASE  ADJ  VN) 


(137) 


(138) 


The  Application  of  R and  K to  Find  Weapon  Radius  - In  the  previous  paragraphs, 
individual  components  of  the  VN  system  are  discussed.  Now  these  components  will  be 


*See  reference  (a). 
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integrated  into  a useful  system.  A target's  vulnerability  may  be  described  by  a VN 
number  consisting  of  3 parts:  a number,  a target  type  symbol  and  a K factor.  The 
number  expresses  the  peak  damage  pressure  (FMO)  for  a given  level  of  damage  for  a 
20-kiloton  weapon,  and  for  a given  probability,  say  .5;  the  target  type  symbol  indicates 
the  most  dominant  damage  pressure  type;  the  K factor  in  effect  adjusts  the  number  (or 
FM)  for  weapon  yields  other  than  20 -kilo tons.  Also,  the  adjustment  factor  as  described 
previously,  is  the  ratio  of  FM  to  FMO  which  can  be  converted  to  an  Adjustment  Value. 
The  Adjustment  Value,  as  explained  previously,  is  applied  to  the  basic  datum  yield  VN. 
Thus,  the  purpose  of  the  whole  K factor  scheme  is  to  derive  this  Adjusted  VN  , winch 
is  the  tool  used  to  determine  the  usable  and  meaningful  value  of  Weapon  Radius  (WR)  for 
a given  yield. 

The  Adjusted  VN  may  be  calculated  by  use  of  the  formulas  just  given.  Having  the 
Adjusted  VN  , one  obtains  the  conveniently  selected  WR^,p  at  the  proper  HOB.  The 

WRiKT  is  then  scaled  to  the  proper  WR  by  means  of  the  cube  root  scaling  law  with  the 

desired  yield. 

WR  = WR1KT  . Y1/3  (139) 


Aiming  Errors,  Target  Distributions,  Kill  Functions  and  the  Determination  of  a Condi- 
tional  Single  Shot  Kill  Probability 

The  preceding  discussion  has  shown  how  the  50  percent  kill  radius  is  determined 
from  semi -theoretical  relationships,  on  the  basis  of  target  hardness  (expressed  as  a 
VN  number  or  in  psi)  and  weapon  yield,  hi  this  subsection,  aiming  errors,  the  distribu- 
tion of  the  target,  and  a set  of  kill  functions  will  be  explored,  and  functions  for  a condi- 
tional single -shot  kill  probability  will  be  derived.  It  is  a conditional  probability,  because 
weapon  delivery  to  the  target  is  assumed.  Removal  of  this  condition  to  find  an  overall 
single -shot  kill  probability  will  be  covered  later. 

Aiming  Error,  Target  Distribution,  and  Impact  Point -Target  Element  Distances  - 
F irst,  consider  the  situation  in  figure  15  in  which  the  targets  are  continuously  distributed 
and  this  distribution  is  circular  normal  with  center  (0,0).  Now,  assume  that  the  center 
of  the  target  distribution  is  also  the  aim  point  for  weapon  delivery  and  that  the  weapon 
delivery  error  is  also  circular  normal  and  centered  about  (0,0).  For  a differential  ele- 
ment of  the  target  area  located  at  (X'^,  X’2)  , the  probability  of  its  being  killed  is  directly 

related  to  its  distance  from  the  differential  element  of  area  at  the  impact  point  (X"j,X"2). 

This  distance  has  components  (XA^,  XA2)  related  to  (X^,  X2)  and  (X'^,  X'2)  as  follows: 


XAl  = X"l  - X^ 
xa2  = x"2  - X’2 


(M0) 


Now,  the  situation  portrayed  in  figure  15  and  equation  (140)  must  be  generalized  if  one 
is  to  find  the  probability  density  distributions  of  XA^  and  XA 2 from  the  probability  density 
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distributions  of  X'^X^,  X’^  and  X"2  . By  the  definition  of  a circular  normal  distribu- 
tion centered  at  (0, 0): 

P(X\)  = l/(cr'  y[2n)expl  - (Xy<7  ’)2/2] 

P(X'2)  = l/(a'  ^2rr)exp[  - (X'1/ct’)2/2 3 
P(X”1)  = l/(a"  ^Tn)exp[  - (X’’1/rf’)2/2] 

P(X”2)=  l/(o"yj2rr)exp[  - (X”2/cr*)2/2  ] 


Impact 

point 


FIG.  15:  AREA  TARGET  SPATIAL  RELATIONSHIPS 


In  equation  (141),  a ' and  a"  are  the  standard  deviations  of  the  distributions.  Gener- 
ally, the  delivery  error  is  specified  in  terms  of  circular  error  probable  (CEP).  This  is 
just  the  radius  of  the  circle  centered  at  the  mean  (in  this  case  (0,0))  which  contains  50 
percent  of  the  impact  points.  For  finding  cr”  from  CEP  it  is  first  necessary  to  find  the 
distance  to  the  impact  point,  R"  . 

R"  = C(X"1)2  + (X”2)2]1/2  (142) 
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Now  that  the  standard  deviations  for  the  equation  (141)  distribution  have  been 
developed,  it  is  possible  to  proceed  to  find  distributions  for  XA^  and  XA2  . It  can  be 

shown  that: 

P(XA1)  = l/i2TT[(a')2  + (a")2]} 1/2  exp  i - (XAl)/' ^2 [(cr*)2  -HO2]}} 

(150) 

P(XA2)  = 1/|2tt[((t’)2  + (a")2])172  exp  j - (XA^Rcr’)2  + (fr")2]}} 

Now,  since  the  probability  of  killing  an  element  of  the  target  is  directly  related  to 
the  distance  between  the  impact  point  and  the  target  element’s  location,  it  would  be  useful 
to  convert  equation  (150)  to  a single  equation  for  the  distribution  of  impact  point  - target 
element  range  R . First:  ■» 


R = [(xa1)2  + (xa2)231/2 


(151) 


As  in  the  previous  cases  of  delivery  error  and  target  distribution,  the  probability 
density  distribution  for  R may  be  written  as  : 

P(R>  = { R/[(<r’)2  + (a’’)2  ]}exp  j - R2/2[(cr’)2  + (a")2  ] J.  (152) 

Substituting  for  a'  and  a"  in  equation  (152),  then  provides  the  desired  distribution 
of  impact  point  - target  element  distances: 

P(R)  = I R/ [(CEP/1, 1774)2  + (RQt./2.4477)2 ] l 

| 2 , 2 2 U <153> 

• exp  j-R7  <2[(CEP/1. 1774/  + (R^/2 . 4477)Z  ] J } 

Generalized  Kill  Functions  and  Weapon  Radius 

A set  of  generalized  probability  of  kill  functions  have  been  postulated  in  reference  (f) 
for  a target  at  a distance  R from  ground  zero.  These  are  of  the  form: 

K“1 

PK(R)  = cxp(  - K • R2/ A2)  I <K  • R2/A2)J/T?  (154) 

J=0 


In  particular,  values  of  K = 3 and  K = 6 in  equation  (154)  have  been  found  to  approxi- 
mate the  dynamic  pressure  and  overpressure  kill  curves  of  reference  (g). 

Since  the  only  parameter  in  equation  (154)  not  defined  is  A , this  will  now  be  derived. 


Weapon  radius  is  defined  in  reference  (g)  as: 
WR  = | y ” 2PK(R)RdllJ-1/2 
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(155) 


X. 


In  equation  (155)  PK  is  the  probability  of  damage  as  a function  of  distance  from 
ground  zero,  R . Now,  by  squaring  equation  (155)  and  substituting  equation  (154)  for 
PK(R): 


(WR)2  =/*  /2R  exp  ( - K • R2/A2)  ^ (K  • R2/A2)  J/J!  1 * 
n * J=0  ’ 


dR 


(156) 


0 

Equation  (156)  can  also  be  written  as: 
K-l  -* 


(WR)2  = £ f 2[(K  • R2/A2)J/J!]exp(  - K • R2/A2)  • R • dR  (157) 

J=0  0 

Consider  the  term  of  equation  (157).  By  making  the  substitution  K • R2/A2  = X, 


rtll 


the  J term  may  be  written  as: 

Jth  TERM  = j'  (A2/KJ!)  * XJ  • exp  ( - X)dX 


(158) 


0 


Carrying  out  the  integration  yields: 
Jth  TERM  = A2/K 


(159) 


Integrating  and  summing  all  K terms  of  equation  (157)  provides  the  following  result: 
9 ^ 9 

(WR)  = E A /K 


J=0 
= A2  . K/K 
= A2 


(160) 


m 


'Plierefore: 

A = WR  (161) 

The  probability  of  killing  a target  at  a distance  R from  ground  zero  is,  then: 

K“1 

PK(R)  = exp  [ - KRZ/(WR)2]  £ [K  • R2/(WR)2  ]J/J!  (162) 

J=0 

K = 3 dynamic  pressure  kill 
K = 6 overpressure  kill 
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The  Probability  of  Kill  for  a Successful  Weapon  Delivery 

The  preceding  two  subsections  have  dealt  with  two  aspects  of  the  kill  problem:  The 
probability  that  a differential  element  of  target  is  a distance  R from  ground  zero  and 
the  probability  of  an  element  of  target  at  a distance  R from  ground  zero  being  killed. 

It  remains  in  this  subsection  to  tie  these  two  together  and  obtain  the  probability  of  kill 
for  a single  weapon  successfully  delivered.  For  this  purpose  conditional  probability  and 
marginal  probability  relationships  will  be  used. 

Consider  two  random  variables  $ and  iji  defined  in  the  same  sample  space  with  fixed 
values  $ = $2  • • • ^ and  ♦ = ♦j*  ^ • The  marginal  probability  of  $ = §j  is 

defined  as  the  probability  that  $ = §j  and  <!)  = ^ summed  over  all  K . Symbolically: 


P<*  = fy  = E P<*  = *T,  ♦ * 

J All  K J R 


(163) 


The  conditional  probability  of  $ = $j  is  defined  as  the  probability  that  § = given 


that  t!f  = . Symbolically: 

P(*  = *j  | * = »K)  = P ($  = , i!f  = <fK)/P(i|t  = *K) 

Combining  both  equations  (163)  and  (164): 

P($  = $ ) = £ P(t  = V)  • P(*  = *T  It  = M 
J All  K K J 1 K 


(164) 


(165) 


In  the  particular  case  under  consideration,  the  probability  of  the  target's  being 
killed  is  equal  to  the  summation  of  conditional  probabilities  of  kill  given  a ground  zero 
distance  of  R times  the  probability  of  ground  zero  being  at  a distance  R . Applying 
tills  to  equations  (153)  and  (162)  and  replacing  the  summation  with  an  integration: 

/*  ® 

PKILLjW  -J  |R/[(CEP/1.1774)2  + (R95/2.4477)2]|  ♦ 

0 

exp  j-R2/|2-  [(CEP/l.  1774)2  + (R95/2.4477)2]U  (166) 

K “1 

| exp  [ - KR2/(WR)2  ] • £ [K.  R2/(WR)2  ]J/J! ! dR 
1 J=0  ’ 


Now,  for  simplicity, 
let  a = K/(WR)2 

b =1/|2 [(CEP/l.  1774)2+  (R95/2.4477)2  ] | 


(167) 
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Equation  (166)  then  becomes: 


PKILL1W  = J «j(2bR  • exp  [ - R2  (a  + b)]  [(aR2)J/J! ]}dR 

0 J=0 


(168) 


rtll 


The  J term  of  equation  (168)  is: 

Jth  TERM  = f (a  Jb/J!)R2J  exp  [ - R2(a  f b)  ] 2RdR 
0 

Integrating  (169)  yields: 

Jth  TERM  = aJb/(a  + b)J+1 


(169 


(170) 


'Hie  complete  series  therefore,  is: 

K-l 

— ru/t* 

'IW 


PKILL  ^ = [b/(a  + b)  ] 5]  [a/(a  + b)  ] J 

J=0 


By  taking  the  sum  of  the  geometric  series  in  equation  (171): 
[b/ (a  + b) ] • 1 1 - [a/(a  + b)]Kj. 


PKILL 


IW 


1 - a/(a  + b) 


= 1 - [a/(a  + b) ] 


K 


(171) 


(172) 


By  multiplying  the  numerator  and  denominator  of  the  term  within  the  brackets  by 
K/ab: 

iK 


PKILL  IW=  1 


r K/2b  T 
L K/2B'+'K72a  "J 


Now,  substituting  equation  (167)  into  equation  (173)  provides: 


PKILL  IW  = 1 - 


K-  [(CEP/1.1774)2  + (R95/2.4477)2] 

K . [(CEP/1.1774)2  + (R95/2.4477)2]+(WR)2/2 


(173) 


(174) 


Of  course,  when  a point  target  is  being  considered,  R^  = 0 in  equation  (174). 
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Strategy  Considerations  and  Weapon  Delivery  Probabilities 

Thus  far,  the  probability  of  kill  development  has  taken  into  account  only  those  factors 
of  significance  for  a weapon  type  actually  delivered  to  the  target.  Now,  the  probability  of 
being  able  to  deliver  a weapon  to  the  target  must  be  considered.  The  overall  single-shot 
kill  probability  of  weapon  type  IW  against  target  type  IT  depends  on  inflight  reliability 
(or  the  probability  that  the  weapon  will  arrive  at  the  area  defenses),  the  probability  of 
penetrating  the  area  defenses,  and  the  probability  that  there  will  be  a target  available  to 
hit.*  Therefore,  the  single -shot  kill  probability  may  be  written  as: 

SSKP  = (INFLIGHT  RELIABILITY)  • (PENETRATION  PROBABILITY)  • 
tvv , n (PR0BAB1L1TY  OF  SURVIVING  ATTACK  BY  RANDOM  TARGET 

TERMINAL  DEFENSE)  • (PROBABILITY  THERE  IS  AN  AVAILABLE 
TARGET)  • (PROBABILITY  OF  KILL  WHEN  WEAPON  IS  DELIVERED) 
(ALLOWABLE  STRATEGY  FACTOR) 


PENPROB. 


PSURViw^t-PTGTjw.it 


PKILL 


STRAT 


IW,  IT 


IW,  IT 
(175) 


The  only  elements  of  equation  (175)  not  previously  considered  are  the  probability 
that  there  is  an  available  target,  the  probability  of  surviving  attack  by  target  terminal 
defense  when  the  terminal  defense  is  random,  and  the  allowable  strategy  factor. 

The  probability  of  surviving  attack  by  the  target  terminal  defense  is  specified  as: 

psurviw,it  = (1  -pkit,iw>  (176) 

The  probability  that  there  is  a target  left  at  the  site  is  the  fraction  of  weapons  with- 
held from  a prior  strike  plus  the  fraction  of  weapons  sent  but  not  successfully  launched. 

The  probability  that  the  target  type  has  not  been  destroyed  previously  is  SURVXjj  . 
Therefore,  the  />robability  that  there  is  a surviving  target  at  the  site  is: 

ptgtit  w = CO  - fsentit)  + fsent]T  • (i  - rrit)  ] survxit 

(177) 

= (1  - FSENTjt  • RRit)  • SURVXit 


Now,  the  strategy  factor,  STRAT^  ^ , enters.  If  the  weapon  type  is  allowed  to 
attack  the  target  type  STRAT ^ ^ = 1 . Therefore,  the  single-shot  kill  probability 
for  a single  warhead  is: 


ssKPijyy,  IT  = rn^ 


PENPROBr 


(1  _PKit,IW^ 


• (1  - FSENTit  • RRit)  • STRATW  IT  • PKILL^  IT 

*For  example,  a bomber  arriving  at  an  ICBM  site  may  find  that  all  ICBMs  have  already 
been  fired. 


Finally,  the  number  of  warheads  per  independently  targetable  weapon  is  used  to  find 
the  single -shot  kill  probability  per  independently  targetable  weapon.  Since  the  number 
of  warheads  (WHDS^)  per  independently  targetable  weapon  (TPRWPN^)  may  be  con- 
sidered as  an  equivalent  number  of  Bernoulli  trials,  then: 

sskpiw,it  = 1 - - sskpiiw  jt)whdsiw/tprwpniw  (179) 

This  completes  the  derivation  of  the  single-shot  kill  probability  relations  hip. 
DAMAGE  CALCULATION  MODELS 

The  damage  calculation  models  are  important,  since  in  many  cases  more  than  a 
single  weapon  will  be  used  against  a target  and  the  single -shot  kill  probability  model 
previously  developed  will,  therefore,  not  be  enough  for  total  damage  calculations. 

There  are  two  different  models  for  calculating  damage.  The  first  is  for  calculating 
damage  to  point  targets,  the  second  for  calculating  damage  to  area  targets  with  a circular 
normal  distribution  of  target  value,  e.g.,  an  approximation  of  fatalities  in  an  attacked 
city. 


Symbols  used  in  the  equations  are  listed  below: 


A 

K 

N 


NP 


IW,1T 


Target  area 

Expected  lethal  area  of  one  weapon 

Number  of  weapons,  or  the  number  of  Bernoulli  trials 

Price  in  independent  targetable  weapons  for  each  weapon  type 


P 

PKILLj  (N) 
R 

RK 

SSKPIW,IT 

survit 

SURVj(N) 


Probability  of  the  weapon’s  being  delivered  to  the  target  area 

Overall  probability  that  the  Ith  increment  of  target  within  area  A 
is  killed  for  N weapons 

Distance  from  the  target  center 

Lethal  radius  of  one  weapon 

Single  -shot  kill  probability  of  weapon  type  1W  against  target  type  IT 
Overall  survival  of  the  target  type 

Overall  probability  of  survival  for  Ith  increment  of  target  within 
area  A for  attack  by  N weapons 


TPKILLjW  j^Total  probability  of  a target's  being  killed 
TVKILL  Total  value  of  target  killed 


VIT 
VKILL 


IW , IT 


Value  given  to  target 
Expected  value  killed  per  target 
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VKILLEDj  Expected  value  of  target  killed  by  N weapons 

W^y  Number  of  weapons  of  type  IW  used  against  target  type  IT 

^*IW  IT  Effective  number  of  weapons 

oo  Weapon  density 

Point -Target -Damage  Model  - No  Terminal  Defense  - Single  Weapon  Type 

In  attempting  to  destroy  a point  target,  a single  successful  weapon  is  assumed  to  be 
all  that  is  required  to  destroy  the  target  completely.  The  probability  of  success  with 
a single  weapon  of  type  IW  against  target  type  IT  is  SSKP^  ^.When  W^  ^ weapons 

of  type  IW  are  used  against  a single  target  of  type  IT  , the  total  probability  of  kill  is: 

w 

tpkilliwit=  1 -(1  -sskpIW(IT)wiw,it  (180) 

Considered  over  many  targets  of  type  IT  the  expected  value  killed  per  target  when 
each  target  has  the  value  Vj,p  is: 

VKILLIW,IT  = VIT[l  ■ 0 -SSKPIWjIT)WIWtIT}  (181) 

This  is  the  final  expression  for  damage  to  point  targets. 

Area -Target -Damage  Model  for  Circular  Normal  Distribution  of  Target  Value  - No 
Terminal  Defense  - Single  Weapon  Type 

For  area  targets  such  as  cities,  where  target  damage  is  measured  in  fatalities  and 
the  lethal  radius  of  a single  weapon  is  small  compared  with  the  size  of  the  target,  the 
point -target -damage  model  is  not  applicable,  since  no  single  weapon  may  be  capable  of 
destroying  the  entire  target  and  if  a single  weapon  is  not  capable  of  destroying  the  entire 
target,  the  weapons  should  be  distributed  in  a way  to  maximize  expected  destruction. 

The  problem  then  becomes  (see  reference  (e)>:  maximize  the  target  value  destroyed 
at  a target  of  type  IT  , subject  to  the  constraint  that  only  W^  ^ weapons  of  type  IW 

may  be  used.  The  solution,  when  obtained,  will  be  in  terms  of  two  limit  curves  for 
expected  damage.  The  lower-value  limit  curve  for  damage  will  then  be  used  for  damage 
calculations  involving  area  targets,  since  it  yields  a more  conservative  result. 

„ 9 

hirst,  assume  that  the  lethal  area,  n(RK)  , of  a single  weapon  is  small  with  respect 
to  the  total  target  area  A . For  a single  weapon,  the  probability  of  the  weapon  being 
delivered  to  the  target  area  A is  P . The  probability  that  an  increment  of  target  at  an 

arbitrary  point  within  A is  killed  by  a weapon  delivered  to  target  area  A is  tt(RK)2/A  . 
The  overall  probability  that  the  Ith  increment  of  target  within  A is  killed  is: 


(182) 


~~  'f  v&z*. 


!>  jJsv?j5/-''v  'v414*=?i'^V 


PKILLjU)  = P • ttRK2/A  (P  • tt  • RK2/A  < 1) 


PKILL^l)  = 1 


(P  • tt  • RK/A  2 1) 


The  corresponding  probability  of  survival  for  the  Ith  increment  is: 
SURVj(l)  = 1 - PKILLj 

= 1 - P • tt  • RK2/A  (P  • tt  • RK2/A  < 1) 


SURVj(l)  = 0 


(P  * tt  * RK/A  a 1) 


For  N weapons  delivered,  corresponding  to  N Bernoulli  trials,  the  probabilities 
of  survival  and  destruction  for  target  increment  I are: 


SURVI(N)=(1  - P • tt  • RK2/A)N 
PKILLj(N)  = 1 - (1  - P . tt  . RK2/A)N 


(P  . tt  . RKVA  < 1) 


SURVj(N)  * 0 
PKILLj(N)  = 1 


(P  . tt  • RK/A  2 1) 


Now,  let  K be  the  expected  lethal  area  of  one  weapon  and  to  be  the  weapon  density. 


K = P • ttRK 


u)=  N/A 

Substituting  equation  (185)  into  equation  (184)  yields: 


SURVjCN)  = (1  - Kco/N)N 


PKILLj(N)  = 1 - (1  - Kto/N)1 


(Kto/N  < 1) 


SURVj(N)  =0 


PKILL  (N)  = 1 


(Kto/N  2 1) 


Now,  assume  that  there  are  M increments  of  the  target,  the  Ith  incremen"  contain- 
ing an  area  AAj  and  having  a value  density  of  Vj  per  unit  area.  The  expected  value 

killed  by  N weapons  is: 
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VKILLEDj  = V • [1  - (1  - Koj/N)in3  aAj  (Ko>/N  < l) 
VKILLEDj  = V AA  (K<o  /N  a i) 


(187) 


For  all  M increments  the  total  value  killed  is: 


M 


TVKILL  = £ V.  Cl“  (1  - Ko)/N)N]AAt  (Ko>/N  < 1) 
1=1  1 1 


M 


TVKILL  = £ VtAAt 
1=1  1 1 


(Koj/N  a 1) 


(188) 


It  is  now  desired  to  maximize  the  total  value  killed  by  adjusting  the  value  of  u‘ 
in  some  manner  corresponding  to  the  value  of  the  target  increments . Thus 

replaces  u)  in  equation  (188).  Now,  introduce  the  constraint  of  total  weapons  expended, 
W , by  the  relationship: 

M 

w = £ 0),  AAt  (189) 

1*1  1 1 

While  it  is  noted  that  W and  N are  actually  the  same,  it  is  useful  in  this  para- 
metric approach  to  keep  them  separate  for  the  theoretical  development. 


The  problem  is  now  to  maximize: 

M vj 

TVKILL  = £ VT[1  - (I  - K<dt/N)  j AA. 

1=1  A 1 1 

M 

TVKILL  = £VTAAT 
1=1  1 1 


(For  KWj/N  < 1) 


(For  KcOj/N  a 1) 


(190) 


Subject  to: 

M 

W = £ a>  AA 
1=1  1 1 

Since  the  first  of  equations  (190)  is  a concave  function  in  N for  each  increment, 
the  total  value  killed  will  be  maximized  when  the  marginal  kills  in  each  increment  are 
equal.  To  establish  the  marginal  kill  for  the  Ith  increment,  both  equations  (190)  are 
differentiated  with  respect  to  . This  yields: 
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d^VKILLVdoJj  = Vj  K(1  - Ku^/Nr"1  AAj  (Ku^N  < 1) 


daVKILLVdojj  = 0 

dW/dcOj  = AAj 


(K  cyN  ;>  1) 


Dividing  the  first  two  equations  of  (191)  by  the  third  provides  the  marginal  kill 
relations  hip: 


d(TVKILL)/dW  = Vj  . K(1  - KWj/N)r 
d(TVKILL)/d\V  = 0 


(Kwj/N  < 1) 


(KoJj/N  a 1) 


Letting  d(TVKILL)/dW  = X , equation  (192)  may  be  solved  for  co. 


x 

wI  = r {'  ' [X/(VI  'K^Wri  } 0/0^  * K)  <1] 


"f0 


[X/Vj  K)S1] 


Substituting  this  result  into  equations  (190): 


M „ 

TVKILL  = E [l  - [X/(^  • K)]1^1  jv^Aj  [For  X/(Vj.  K)  < 1 ] 


TVKILL  = ^ V.  A A. 
1=1  1 


w-z  £{i-  [x/(v:  . k)]1^}  aAj 


[For  X/(Vj  • K)*l]  (194) 


[For  X/(V j ♦ K)il] 


Now,  for  area  targets,  the  value  is  assumed  to  have  a circular  normal  distribution, 
in  which  R is  the  distance  from  target  center  and  is  the  standard  deviation: 

V = [l/(2no,2y)]  exp  [ - R2/(2cr2v)]  (195) 

2 

Rewriting  equation  (195)  and  expressing  R in  terms  of  V yields: 


R2  = -2(?2v  loge(2TTayV) 


V s l/(2nff  y) 
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Multiplying  equation  (196)  by  n would  yield  the  enclosed  area  as  a function  of  the 
enclosed  value.  By  differentiating  the  results  dA  may  be  obtained  in  terms  of  V : 


d(nR2)  = dA  = - 2ttct2v  dV/V 


(197) 


Now,  allowing  the  area  increment  to  decrease  to  differential  size  and  allowing  M-*», 
equation  (194)  may  be  written  as: 

N 

TVKILL  |l  - [X/(V  • Kjf1^  } VdA 


W =/  K { 1 " [X/V  * *>3^}  dA 


J 

1 


(198) 


Substituting  equation  (197)  into  equation  (198)  and  using  the  limits  dictated  by 
equation  (194)  and  (196)  provides: 


1/(2ttct  y) 


N 


TVKILL  = 2tt  a 


v/ 

X/K 


|l  - [X/tV.KJl^jdV 


J 


W = 2tt  a 


l/(2nor  y) 

2f 


V 
X/K 


£ [l  - O/V  • K)]^-^ 


Integrating  equation  (199)  yields: 


N 


TVKILL  = 1 - 2tto r2y  X/K  + (N-l)  [^na^X/K)*^  - 2n<y2vX/K] 

W = (2na2vN/K)  j-  loge(2TT<72vX/K)  + ^na^X/K)^-!  ](N-1)| 


(199) 


(200) 


Now,  since  N and  W have  been  kept  separate,  it  is  possible  to  define  two  limit 
curves  for  the  extreme  conditions  of  N,  namely  N-*l  and  N-»».  These  two  conditions 
imply  the  situation  of  a single  weapon  or  an  infinite  number  of  weapons  assigned  to  the 
target.  To  find  the  equation  applicable  to  these  two  limiting  cases,  however,  requires 
some  effort  since  by  inspection  it  can  be  verified  that  both  of  equation  (200)  have  inde- 
terminate forms  for  N=1  and  N = » . 

First,  in  equation  (200)  the  substitution  is  made: 
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3 = (2m/y\/K)  1N_i 

Equation  (200)  may  now  be  rewritten  as: 


TVKILL  = 1 - p -i[l  + (N-l)(l-/B)] 


2tkt  yN(N-l) 


C3  - 1 - logg(3)  ] 


The  second  of  equation  (202)  may  also  be  rearranged  to  yield: 


2ttct  yN(N-l) 


= 3 ~ 1 - log  0)  = X 


By  noting  the  original  limiting  conditions  in  equation  (199)  and  elsewhere,  it  may  be 
declared  that  the  range  for  3 is  (0  sjl  s 1).  For  the  limiting  values  of  3=0  and  3=1, 
it  is  useful  to  tabulate  the  corresponding  values  of  N and  X from  equations  (201),  (202) 
and  (203): 

3 N X 

0 1 ® 

1 “ 0 

Thus,  the  conditions  of  3=0  and  3=1  correspond  to  the  condition  of  N-*l  and 
N-»»  , respectively,  the  conditions  of  interest  as  stated  earlier. 

The  N -*1  case  will  now  be  listed  in  detail. 

Rewriting  equation  (203),  an  expression  for  N-l  may  be  obtained: 

N-l  = C/(NX) 

9 (204) 

C = KW/(2tt(t"v)  ' ; 


Using  equation  (204)  in  the  first  of  equation  (202)  results  in: 

C 

TVKILL(l)  = 1 -3"^{l+  [C/(NX)][1-3]| 

C 

= 1 - 3 NX  |l  + C/(NX)  - 3C/(NX)  j> 


Now,  as  established  previously,  as  N-*l,  P -»0  and  X-*“.  Therefore  the  last 
two  terms  go  to  zero  as  N -»1  and: 


TVKILL(l)-»l~j3 


(206) 


Equation  (206)  may  now  be  rewritten  and  the  natural  logarithms  of  both  sides  taken: 


loge[l  - TVKILL(l)]  = (C/X)  loge((3) 


(207) 


Using  the  equivalent  expression  for  X from  equation  (203)  in  equation  (207)  results 


loge[l  -TVKILL(l)]  = Cloge(^)/[j3-i  -log^)] 

Now,  in  the  limit  as  N -*1 , (P  ->  0): 

loge[l  - TVKILL(l) ] -»c/l/[  - 1 + 0-l)/loge0)]| 
(3-.0 

-♦  -C 


(208) 


(209) 


The  n: 


1 - TVKILL(l)  = exp(-C) 


= exp(  - KW/(2ttct  v)) 


(210) 


Therefore,  the  N=1  bound  of  the  circular  norma]  target  destruction  easels: 

TVKILL(l)  = 1 - exp(  - KW/(2na2y)) 

The  N -♦  “ case  will  now  be  treated  in  detail. 

First  of  all,  when  N -»  ®>,  equation  (203)  maybe  approximated  by: 


2ncr  yN(N-l)  2na  yN 


N = (C/X)J 


C = VW/(2na  v) 
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By  substituting  equation  (213)  in  the  first  of  equation  (202): 


TVKILL(®)  = 


1/2  -1 

jl  + [(C/X)1/2 


l](H3)j 


= 1 - /3 


1/2 


[1  +(C/X)1/2(1  - 1/0)] 


(214) 


The  last  term  in  equation  (203)  is  indeterminate  as  N -»»  , since  X-»0  and  0-»l  . 
The  limit  for  this  last  term  may  be  found,  however,  by  substitution  of  the  equivalent 
expression  for  X from  equation  (203)  and  then,  for  ease  of  operation,  evaluation  of 
the  limit  of  the  last  term  squared  as  N -*<*>  (0  -»1). 


LAST  TERM  = (C/X)1/2(l  - 1/0) 

_ C1/2(l  - 1/0) 

[0-1  - loge(0)  ]1/2 


(215) 


Then: 

(LAST  TERM)2  = g'.1 


(216) 


Differentiating  numerator  and  denominator  witn  respect  to  0 (1  'Hospital's  Rule) 
to  find  the  limit  of  equation  (216)  as  0-*l  yields: 

(LAST  TERM)2  -»  2C(-  " (217) 

0-*l  1-1/P 

-♦  2C 

Therefore: 

LAST  TERM  -»  V2C)1/2  (218) 

In  equation  (207)  the  positive  root  was  selected  since  C , (1-0)  and  X are  all  ^0  . 
Now  equation  (203)  has  become: 


TVKILL(»)  = 1 - 


1/2 

[1  + (2C) 


(219) 
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Comparative  kill  results  from  equation  (211)  for  N-*l  and  from  equation  (224)  for 
N -4  «are  shown  in  figure  16.  Since  the  kill  values  obtained  by  using  TVKILL(«)  are 
always  less  than  the  values  obtained  by  using  TVKILL(l),  TVKILL(»)  will  provide  a 
conservative  estimate  of  kill  and  has  been  used  as  the  basic  model  for  evaluating  kill 
for  area  targets. 


Although  equation  (224)  has  incorporated  weapon  delivery  probability,  delivery  error 
has  not  been  included.  It  therefore  remains  to  integrate  this  into  the  damage  model.  To 
do  this,  the  value  of  W is  set  equal  to  one  in  equation  (224)  and  the  result  set  equal  to 

the  single-shot  kill  probability  with  K'  used  in  place  of  [K/(2rr<72.y)]^2.  Thus: 


SSKP1W,IT  = (1  + K'lW,lT)  exp  ( ' K IW,IT)  (225) 

The  resulting  value  of  K'  is  then  found  by  an  appropriate  method  (such  as  Newton's 
method).  The  final  damage  model  for  area  targets  is  therefore: 

VKILLj^  = VIT[l-tK’IW  1T(Wj^  IT)  ]exp[-  K'IW> iX^WIW, IT^  / 1(226) 


Using  this  approach  the  damage  from  equation  (226)  and  from  equation  (181)  as  a 
function  of  the  number  of  weapons  employed  would  be  as  shown  in  figure  17.  It  can  also 

2 

be  shown  that  for  some  ^ (in  place  of  K/(2rrcr  v))  determined  from  SSKPjW  ^ , 

equation  (181)  and  equation  (211)  yield  identical  results  for  the  same  single  shot  kill 
probability.  Thus,  the  optimistic  area  target  damage  model  and  point  target  damage 
model  are  identical  for  identical  single-shot  kill  probabilities. 


Point-Target-Damage  Model  - Terminal  Defense  - Multiple  Weapon  Types 

In  this  extension  of  the  point -target -damage  model  developed  previously,  it  can  be 
considered  that  each  weapon  type  must  pay  a price  at  the  target  of  NP^y  ^ independently 

targetable  weapons.  For  Wp^  ^ weapons  of  type  IW  allocated  to  target  type  IT, 

the  effective  number  of  weapons  is: 

c 


W*  = W - NP 

IW,  IT  IW,IT  IW,IT 


(227) 


The  overall  survival  of  the  target  type  is,  then,  the  probability  that  it  survives  each 
weapon  type’s  raid  or: 


SURVtt=  n SURVrr 

IT  ALL  IT,IW 

WEAPON 

TYPES 

ASSIGNED 


(228) 


Since: 

pkillit  = i - survit 


(229) 
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Then  equation  (181)  may  be  extended  by  using  equations  (227),  (228)  and  (229): 

vkill1t»vit[i  - n (i-sskpIWjIT)wiw,it“npiw,it 

weapon 

TYPES 

ASSIGNED 


(230) 


Area -Target -Damage  Model  - Terminal  Defense  - Multiple  Weapon-Types 

In  extending  the  area -target-damage  model,  the  process  is  not  so  direct  as  in  the 
point -target  case.  While  it  is  true  that  equation  (227)  applies  in  arriving  at  the  effective 
number  of  weapons,  equation  (228)  does  not.  Presumably,  the  total  number  of  weapons 
should  be  applied  in  some  optimum  density  pattern  over  the  target.  How  the  parameters 
can  be  combined  can  best  be  seen  by  consideration  of  two  weapon  types  with  identical 
kill  probabilities  against  target  type  IT  . Then  W . weapons  of  the  first  type  and 
^2  IT  weaPons  a second  type,  in  the  absence  of  a terminal  defense  should  yield  the 

same  results  as  if  jy  + W£  jy  weapons  of  either  type  had  been  used.  This  leads 

by  heuristic  argument  to: 

VKILLIT  = VIT|l +[^(K'IW  IT)  ^WtlT 

^ ALL 
WEAPON 
TYPES 
ASSIGNED 

• {-p{-CBk'iw,ix)2w1w 

V.  A I I 

WEAPON 
TYPES 
ASSIGNED 

When  the  terminal  defense  is  taken  into  account  by  use  of  the  effective  number  of  weapons 
W*jw  jy  from  equation  (227),  equation  (231)  becomes: 

VKILLjy  = VjT  jl  +[£(K'I\y,rr)  (WIW,IT  ” NPIW,IT^  / f 

t.  ALL  J 

WEAPON 
TYPES 

ASSIGNED  -s  (232) 

• {expfCBK'iw.iT)  <Wiw,IT  ' NPIW,IT^  1 }} 

WEAPON 

TYPES 

ASSIGNED 


This  completes  the  derivation  of  the  damage  calculation  models. 


Sgg^ggggggggaiaBffjBj 
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WEAPON-TARGET  ALLOCATION  MODEL 

The  efficient  allocation  of  resources  has  always  been  one  of  the  core  problems  of 
economic  theory.  To  solve  such  problems,  a great  body  of  theoretical  literature  in  the 
fields  of  linear  and  non-linear  programming  has  been  generated  in  recent  years.  One 
set  of  methods  within  this  body  of  literature  makes  use  of  the  so-called  LaGrange 
Multiplier  Theory  as  a convenient  means  of  solution  for  resource  allocation  problems. 

The  basic  resource  allocation  problem  to  be  solved  in  Code  50  by  the  weapon-target 
allocation  model  is  the  allocation  of  weapons  (resources)  to  targets  so  as  to  maximize 
the  destruction  of  enemy  target  value  subject  to  the  constraints  of  weapon  inventory, 
integer  number  of  weapons  assigned  to  each  target,  and  lower  limits  for  the  destruction 
of  defended  targets  (if  achievable).  The  solution  is  provided  by  a procedure  making  use 
of  LaGrange  Multiplier  Theory,  steepest  ascent  methods  of  mathematical  programming 
and  certain  heuristic  mathematical  programming  algorithms  which  have  proved  sufficient 
in  obtaining  a solution  to  this  problem. 

Symbols  used  in  the  equations  of  this  section  are  listed  below: 


CRIT 

critit 

CRITFAC 

°IT 

DPREF 

DSUM 

DTOT1 

dtot2 

FRACINV 

FRACP 

F RAC PA 
FRL 

FRU 

HIW,IT 


Maximum  value  for  CRIT^ 

Penalty  parameter  for  target  type  IT  for  using  the  second  best 
weapon  type 

Assured  destruction  modifier  for  CRIT  based  on  whether  assured 
destruction  criteria  have  been  satisfied 

Damage  to  target  type  IT  due  to  weapon  types  W^  ^ WN 

assigned  to  target  type  IT 

Total  defended  target  value  destroyed 

Total  damage 

Initial  allocation  of  damage 
Second  allocation  of  damage 
Class  of  weapons 

Fraction  of  defended  target  value  specified  for  destruction;  also 
called  the  preferred  target  destruction  criterion 

Fraction  of  defended  target  value  destroyed 

Fraction  of  defended  target  value  destroyed  for  undefended  target 
multiplier  value  of 

Fraction  of  defended  target  value  destroyed  for  undefended  target 
multiplier  value  of 

Number  of  weapons  which  maximizes  the  LaGrangian  profit  for 
single  target  type 


-71- 


HMAXCjT 

HMAXNit 

HSUM 

ITC 

IWM 

KIW,IT 

MAXWPS 


NIIWfIT 

NOTAS 

NPIW 
NWP, 


IW 


IW 
NWPD 


rw,  it 


PERR 


I 


PERR, 


PERR, 


sskpiw,it 

STOGO 

SVAL 

SVDL 

UBOUND 

VIT 

VPREF 


W* 

IW,  IT 
W* 

w IWM,  ITC 


Upper  bound  of  the  LaGrangian  profit  for  target  type  IT 

Best  LaGrangian  profit  per  target  for  target  type  IT 

Second  best  LaGrangian  profit  for  target  type  IT 

Overall  upper  bound  to  the  LaGrangian  profit 
Most  critical  class  of  targets 
Best  type  of  missile 

Expected  lethal  area  of  weapon  type  IW  used  on  target  type  IT 

Number  of  weapons  applied  which  are  of  best  type  for  use  against  at 
least  a single  target  in  the  critical  class 

Integer  number  of  weapons  for  each  weapon  type  that  will  maximize 
the  expected  value  destroyed  per  weapon 

Number  of  unassigned  weapons 

Price  in  terms  of  independent  targetable  weapons  successfully  launched 

Total  number  of  independent  targetable  weapons  available 

Total  desired  inventory  for  weapon  type  IW 

Percentage  e’ror  between  initial  allocation  of  damage,  DTOT^  , 
and  upper  bound  of  damage,  UBOUND 

Percentage  error  between  second  allocation  of  damage,  DTOT^  , 
and  upper  bound  of  damage,  UBOUND 

Percentage  error  between  second  allocation  of  damage,  DTOT2, 
and  upper  bound  of  damage,  UBOUND,  based  on  revised  shadow  values 
Single  shot  kill  probability  for  target  type  IT 

Total  shadow  value  of  the  unassigned  weapons 

Shadow  value  of  the  weapon  inventory 

Shadow  value  for  the  aggregate  desired  inventory 

Upper  bound  for  target  damage  considered  over  the  target  set 

Value  of  target  type  IT 

Total  defended  target  value 

Optimum  number  of  weapons 

Optimum  number  of  weapons  to  be  used  against  a single  target  of  the 
critical  type 
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XMAXL  Absolute  value  of  the  maximum  adjustment  of  X ^ 

a ^ Highest  pr:,>r  value  of  to  produce  an  acceptable  level  of 

assured  destruction  ' 


a , 


N 


a, 


U 


AX' 
ANWPD 

XIWM 


Undefended  target  value  multiplier 

Lowest  value  of  to  produce  an  unacceptable  level  of  assured 
destruction  1 

Correction  factor 

Inventory  increment  for  weapon  type  IW  against  target  type  IT 
Shadow  value  for  the  best  weapon 


LaGrange  Multipliers  from  the  Classical  View 

The  mathematical  formulation  of  the  allocation  problem  may  be  stated  as: 

Maximize: 


y = yOvv-'-V 

subject  to  the  constraints: 

^1^X1*  X2’  * * ,XN^  = ® 
^^i*  x2*  • * *xj^  = ^ 

W'V  X2’  • • *x[vj)  = 0 


(233) 


It  has  been  shown  in  references  (h)  and  (i)  tliat  when  . . .x^  are  continuous 

variables,  y^^, . . .x^)  is  a concave  function  and  M s N (there  are  no  non -binding 

constraints)  a global  solution  to  the  above  problem  can  be  found  by  finding  a solution 
to  the  unconstrained  problem: 


Y = y(xrx2> 


.x 


N 


M 

)-  £ 

K=1 


xk*k(x 


l’X2’  ” 


•*N> 


(234) 


where  the  X^'s  are  non-negative  multipliers.  The  solution  to  equation  (234)  is  then 
generated  by  differentiating  Y with  respect  to  each  and  setting  the  derivatives  to 
zero.  Thus  solving  the  M constraint  equations  and  the  N equations 


M 

3Y/AXj  = dy(x^, X2»  ...x^j)/9Xj  - Xj^(3$j^/dXj)  — 0 

K—  1 

I = 1,2,  ...N 


(235) 
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results  in  the  desired  solution,  a set  of  x°j"s  and  X^’s  which  yield  a y(x°^,x°2, ..  •xV 
that  is  a maximum. 

As  Samuelson  noted  in  reference  (j),  the  values  of  the  X^’s  in  equation  (234)  are 

equal  to  marginal  costs  or  a unit  price  for  each  of  the  M facilities.  However,  since 
these  are  not  necessarily  real  prices  in  problems  of  tins  type,  the  X^’s  are  considered 

’’shadow  prices".  The  interpretation  of  the  summation  on  the  right  side  of  equation 
(234)  is  then  the  "shadow  value"  expended.  The  composite  view  of  equation  (234),  then, 
is  that  the  difference  between  the  "payoff*  and  the  "shadow  value"  expended  may  be 
termed  a "profit".  It  is  this  "profit"  that  is  to  be  maximized. 

The  Generalized  LaGrange  Multiplier  Theory 

For  the  allocation  of  weapons  in  Code  50,  the  classical  problem  may  be  stated  as: 

Maximize  the  total  damage  DSUM: 

DSUM  =£DIT(W1>  it.  . .WN>  JT)  (236) 

ALL 

TARGETS 

IN 

TARGET 

SET 

Subject  to  the  N constraints: 

ZWjw  it  sNWPrW>  IW  = 1,2,...N 

TARGETS  (237) 

IN 

TARGET 

SET 

WIW,  IT  = illteSers 

and  subject  to  the  strategy  constraints  of  no  mixing  of  weapon  types  on  terminally  defended 
targets  and  a minimum  level  of  damage  (if  achievable)  of  x percent  against  that  subset 
of  targets  with  terminal  defenses . 

The  classical  theory  would  then  indicate  that  an  augmented  equation  of  the  form  of 
equation  (234)  would  be  formed.  Then  the  augmenteu  equation  would  be  differentiated 
with  respect  to  each  jT  and  the  result  solved  for  the  weapon  allocation.  This 

procedure,  however,  is  not  tractable  in  this  case  because  (1)  the  numbers  of  weapons 
allocated  to  individual  targets  must  be  integers  and  (2)  the  form  of  the  destruction 
equations  from  the  previous  subsection  is  such  that  solutions  are  not  directly  available 
by  these  means.  To  get  around  this  difficulty,  Everett,  in  reference  (k),  showed  that  the 
theory  of  LaGrange  Multipliers  is  not  limited  in  usefulness  tq  differentiable  functions. 
Essentially  Everett's  method  is  based  on  two  theorems.  These  two  theorems  are  related 
to  the  use  of  the  LaGrange  Multipliers  and  to  the  computation  of  an  upper  bound  with 
which  to  test  a particular  allocation.  These  theorems  are: 
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Theorem  1:  Let  S be  the  domain  of  an  objective  function  F(x)  . If  X^,  K=1,...,N 
are  non-negative  real  numbers  and  if  x*  e S maximizes  the  function: 

K*N 

H = F(x)  - £ X„C„(x)  (238) 

K=1  K K 


Where  the  Cj/s  denote  the  constraints,  overall  x eS,  then  x*  maximizes  F(x)  over 
all  those  x eS  such  that  C^(x)  s C|,(x*)  for  all  K,  regardless  of  the  nature  of  the  set  S . 

As  demonstrated  earlier,  a LaGrange  Multiplier  problem  is  normally  solved  by 
differentiating  the  LaGrangian  function,  setting  the  results  equal  to  zero,  and  solving 
for  both  the  X^'s  and  the  resource  allocation  that  maximizes  the  objective  function, 
given  a set  of  constraints.  This  theorem  says  that  if  the  LaGrangian  is  maximized, 
given  a set  of  X^'s  (i.e.,  given  the  shadow  prices  of  the  resources),  the  resulting 

allocation,  x*  , implies  a set  of  constraints,  C^(x*)  . Thus,  automatically  there  is 

an  optimum  allocation,  x*  , for  a constrained  problem,  namely  the  problem  with 
constraints  CK(x*)  . 


Thus,  there  is  a great  advantage  in  starting  with  the  multipliers  and  deducing  the 
constraints,  rather  than  deducing  the  multipliers  given  the  constraints,  in  that  this 
approach  allows  one  to  place  no  restriction  on  the  domain  of  the  objective  function. 
Therefore,  having  an  objective  function  that  is  neither  continuous  nor  differentiable  no 
longer  presents  an  insurmountable  problem. 

The  disadvantage  of  this  modified  method  is  that  unless  the  correct  set  of  X^’s  is 

chosen  all  the  constraints  will  not  be  satisfied.  However,  one  can  prove  that  if  all  but 
one  constraint  are  held  constant,  the  remaining  constraint  is  a monotone  decreasing 
function  of  its  associated  multiplier.  Thus,  after  a wrong  solution  has  been  reached, 
the  direction  to  adjust  the  X^'s  in  order  to  obtain  the  desired  constraints  is  known. 
Algorithms  have  been  developed  which  make  this  adjustment  relatively  efficient. 

In  many  cases,  one  cannot  find  a set  of  multipliers  which  imply  the  desired  con- 
straints. The  above  theorem  insures  only  that  if  one  finds  a solution,  it  is  a correct 
solution;  it  does  not  guarantee  that  a solution  can  be  found.  The  existence  of  optimum 
solutions  that  can  be  found  by  adjusting  the  multipliers  depends  upon  the  objective  function 
being  approximately  concave  in  the  region  of  the  solution.  Such  concavity  is  not  always 
present;  however,  the  following  theorem  allows  bounds  to  be  placed  on  the  solution  to 
problems  where  a direct  adjustment  of  the  multipliers  will  not  work. 


Theorem  2:  If  x'  comes  within  e of  maximizing  the  LaGrangian,  i.e.,  ifforallxeS, 
where  S is  the  set  of  all  possible  strategies 


F(x')  - £XkCk(x')  > F(x)  -I>kCk<x>  - e , 


(239) 


then  F(x')  , the  payoff  at  x'  , is  within  e of  the  maximum  payoff  for  those  constraints. 
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The  proof  of  this  theorem  gives  some  insight  into  its  usefulness . Rearranging  the 
above  inequality: 

F(x)  - F(x’)  +V;\k[Ck(x')  - CK(x)  3 < e (240) 

Keeping  the  X^’s  fixed,  the  allocation  x*  implies  constraints  C^(x')  and,  obviously, 
x'  is  a feasible  solution  to  the  problem: 

Maximize:  F(x) 

subject  to: 

CK(x)sCK(x')  (241) 

For  any  other  feasible  solution  to  this  problem,  C^(x)  <;  C^(x’)  by  definition,  implying 

CK(x')  - C^(x)  a 0 . Thus,  substituting  into  inequality  equation  (239)  above, 

F(x)  - F(x’)  < e for  any  x that  is  a feasible  solution  to  the  problem  defined  by  equation 
(241)  above.  In  practice,  the  value  of  e is  restricted  by  imposing  the  constraint  that 
the  number  of  weapons  available  for  allocating  to  any  one  target  is  restricted  to  the 
number  of  weapons  in  the  stockpile.  No  payoff  can  be  more  than  e greater  than  the 
payoff  given  by  the  allocation  x'  . 

Thus,  this  theorem  allows  an  upper  bound  to  the  maximum  payoff  to  be  computed 
for  any  constraints.  A set  of  X^'s  is  picked  and  both  the  maximum  LaGrangian  and 

the  implied  constraints  are  computed.  A new  allocation,  x ' , is  then  found  that  implies 
the  desired  constraints,  and,  keeping  the  X^'s  fixed,  the  LaGrangian  at  x'  is 

evaluated  e is  then  the  difference  between  the  maximum  LaGrangian  and  the  value  of 
the  LaGrangian  at  x = x'  . This  is  the  maximum  additional  payoff  above  F(x’)  which 
could  ever  be  obr  aed  while  meeting  the  desired  constraints.  Geometrically,  Theorem  2 
says  that  the  hyperplane  with  slopes  defined  by  the  X^'s  , tangent  at  xq  to  the  envelope 

of  the  maximum  payoff  surface  in  the  space  of  maximum  payoff  versus  resources  ex- 
pended, is  an  upper  bound  to  F(x)  at  any  point  if  xq  maximizes  F(x)  -^X^Cj^x)  . 

The  Mathematical  Programming  Procedure  in  Code  50* 

The  mathematical  programming  procedure  in  Code  50  for  weapon  allocation  involves 
three  fairly  separate  steps.  These  are:  (1)  the  initial  allocation,  (2)  testing  of  the  initial 
allocation,  and  (3)  computation  of  a second  and,  if  necessary,  a third  allocation  if  the 
prior  allocation  is  determined  to  be  unsatisfactory.  These  steps  will  be  discussed  in 
turn. 


*The  remainder  of  this  section  parallels  the  discussion  in  reference  (1). 
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Tlie  Initial  Allocation  - The  first  of  these  allocations  is  a heuristic  process  malting  use 
of  a steepest  ascent  algorithm . The  various  operations  involved  are  listed  and  discussed 
in  the  order  in  which  they  are  carried  out. 

The  first  operation  that  is  carried  out  is  to  find  the  integer  number  of  v/eapons 
NIjw  u for  each  weapon  type  that  will  maximize  the  expected  value  destroyed  per 

weapon.  Of  course,  for  targets  without  terminal  defenses  a single  weapon  produces  this 
result  since  the  destruction  functions  (see  figure  16)  are  concave  functions  and  each 
additional  weapon  after  the  first  adds  progressively  less  value  destroyed.  For  defended 
targets  the  number  of  weapons  producing  this  result  is  some  number  greater  than  the 
price  paid  to  the  terminal  interceptors.  When  the  price  to  be  paid  to  terminal  interceptors 
exceeds  the  number  of  weapons  in  inventory,  the  target  cannot  be  attacked  by  that  weapon 
type.  These  relationships  are  summarized  in  equation  (242). 

NIIW,IT  = 1 : NPIW,IT  = 0 

niiw,it  3 dit(niiw,it)/nim,it  isMAXIMIZED  (242) 

NIIW,rr  •=  • :NPIW>rr>NWPIW 


Now  the  best  payoff  target  class  for  the  first  weapon  type  is  selected.  And  a maxi- 
mum of  1/5  of  the  weapons  of  this  type  are  assigned  to  the  target  class  provided  this  is 
enough  to  attack  at  least  one  target  of  the  class  under  the  maximum  efficiency  attack 
relations  hips  established  in  the  first  two  of  equation  (242).  If  1/5  of  the  weapons  cannot 
produce  a maximum  efficiency  attack  on  a single  target  of  the  class,  then  enough  weapons 
are  allocated  to  the  class  for  attacking  one  target.  Destruction  for  the  target  class 
attacked  is  then  computed.  When  all  targets  of  the  class  are  net  attacked,  the  class  is 
split  into  subclasses  — the  targets  of  the  class  not  attacked  with  the  weapon  type  forming 
the  new  subclass.  This  procedure  continues,  all  target  subclasses  being  considered  for 
remaining  weapon  increments  (1/5  inventory  or  number  for  maximum  efficiency  attack 
on  one  target  of  highesc  payoff  target  subclass  remaining  as  outlined  above)  of  the  first 
weapon  type  until  all  weapons  of  tliis  type  have  been  allocated.  The  procedure  above  is 
then  continued  (except  terminally  defended  subclasses  already  attacked  with  another 
weapon  type  are  not  considered)  for  the  second  weapon  type,  l'-°n  the  third  until  all 
types  have  been  allocated. 


Testing  of  the  Initial  Allocation  - After  all  weapons  of  all  types  iiave  been  assigned  the 
maximum  additional  payoff  available  from  an  addition,  t weapon  of  eacli  type  is  computed, 
witli  the  stipulation  that  no  weapon  mixing  against  terminally  defended  targets  is  allowed. 
This  maximum  incremental  payoff  can  be  considered  the  value  of  a single  weapon  in  its 
best  alternative  use  and  in  accordance  with  economic  theory,  this  is  its  value,  or  as  used 
here  the  "shadow  value"  of  the  weapon  type.  In  mathematical  notation  using  a difference 
eiuation  approach  (since  integer  numbei  of  weapons  are  considered)  and  noting  that  the 
"shadow  value”  is  the  current  value  of  the  LaGrange  Multiplier,  X, 


■iw 


X 


m = MAX<ADSUM>  = 


MAX  <D 

OVER  ALL 
PERMISSIBLE 
tARGLr  CLASSES 


(W 


IW)IT+  1)  - dit(wIWjIT)> 


(243) 
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This  result  also  approximates  the  result  which  might  have  been  obtained  by  formally 
differentiating  the  augmented  destruction  equation,  of  the  form  of  equation  (234)  with 
respect  to  ^ . Noting  that  the  augmented  function,  exclusive  of  strategy  con- 

straints is  of  the  form: 


Y ^ °IT(W  1 , IT'  W2,  IT*  * *' WN,  IT* 


W^wiw,it  “ m?m 


(244) 


differentiating  equation  (244)  with  respect  to  each  W ^ and  setting  the  derivative 
e*  lal  to  zero  results  in: 


d[DIT(W^IT,W2>IT,  •••wNjIT)l/dWIW>IT  - = 0 


(245) 


If  the  derivative  in  equation  (245)  is  approximated  by  a difference  in  which  dW^y 

= 1 , then  it  is  seen  that  the  two  results  evaluated  for  this  allocation  are  approximately 
equal. 


The  formal  evaluation  process  now  continues  by  evaluating  how  much  LaGrangian 
profit  (damage  minus  weapon  shadow  value)  would  result  from  employing  the  best  weapon 
type  on  each  target  (provided  that  the  total  inventory  would  not  be  exceeded  for  any  one 
target  of  the  class).  This  provides  an  "upper  bound"  to  the  LaGrangian  profit  against 
which  the  initial  allocation  may  be  evaluated.  If  the  upper  bound  and  the  initial  allocation 
fall  within  a specified  amount  of  each  other,  the  allocation  process  can  be  terminated  at 
this  point  and  the  initial  allocation  adopted  as  the  final  allocation.  The  upper  bound 
determination  is  based  upon  pure  weapon  allocations,  i.e.,  a single  weapon  type  per 
target  of  each  type.  The  computation  of  the  upper  bound  is  based  upon  the  steps  outlined 
below. 


For  each  weapon  type,  the  number  of  weapons  that  maximize  the  LaGrangian  profit 
Hit  IW  *s  ^ounc*  ^or  a s*ngle  target  of  that  type: 

hit,  iw  = Drr*wiw,  it*  " x iwwiw,  it  (246) 

For  point  targets  and  area  targets  Dj^fW^y  ^)  is: 

Point  targets: 

dit(wiw  , it}  = vit{  1 ■ (1 " SSKPIW , Wrv  ’ IT  " N?IW  ’ IT  } 

Area  targets:  (247) 

^T^IW.IT^  = VIT  j1  “ {X  + ^WIW,IT  " NPIW,IT^KIW,IT*  ]1/2} 


exp]  - ■ npiw,it^kiw,it)  ^ / 1 
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»"  “tins  " (246>  tor  - "»  •"»•  * -geB  and  sottl„g  tha  results  to 
Point  targets: 

dHT 


IT 

m v 


IW.IT 


IT{tloge(l  -SSKPIW|IT)](1  -SSKPIWiIT,WIw  n.-NPiw>iT  1 


Area  targets: 
dH 


" XIW  = 0 


IT 


' =XP  -{C(KIW> 


(248) 


■T>  wiw,it-^w,,t»11/2} 


XIW  =0 


Using  the  form  of  equation  f94fU  n ....„ 

previously  found  i„  equation  (243)  the  . M “ ‘l,e  tyi’e  a"d  the  value  ot  xm 

found  as:  <243>'  t,K  ■“*«  0'  weapons  w . may 


Point  targets: 


W * 
IW.IT 


= I"  W^viTlog,a  - sskpiw  TT)]j 


Area  targets: 


ww,it*  = npiw,it  + 


The  value  of  W 


IW,  IT' 


f lQge  ^VIT(KIW . IT)2/ <2  * nv)  2 

i — p — — ill 

1 IW.IT 


(249) 


(247). 


IW ,IT*  can  then  be  used  to  find  H * frnm 

IT,  IW  from  equations  (246)  and 


HIT,IW*  " “iT^IW.IT*)  • AIWWIW,it' 


Th  (250) 

Value0f  HIT,  IW*  Produced  by^any  sin^^  Therefaref  IT  n,«xlmum 
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OVER 

ALL 
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CLASSES 


IT,IW 


*> 


(251) 


Extending  this  to  the  complete  target  set,  the  overall  upper  bound  to  the  LaGrangian 
profit,  HSUM,  is  just  the  sum  of  all  such  H^*  ’s  . 


HSUM  = E H.  * 

IT 


(252) 


Now  HSUM  is  the  difference  between  the  upper  bound  of  target  damage  UBOUND 
considered  over  the  target  set  and  the  shadow  value  SVAL  of  the  weapon  inventory. 


HSUM  = UBOUND  - SVAL 


SVAL  = E X 
IW 


NWP 

iw1NWtiw 


(253) 


Therefore 

UBOUND  = HSUM  + SVAL  (254) 

The  expected  total  damage  to  the  target  set  provided  by  the  initial  allocation, 

DTOTi,  may  now  be  compared  with  the  upper  bound  of  damage  UBOUND,  using 
Theorem  2 and  the  difference  e computed. 

Gj=  UBOUND  - DTOTx  (255) 

This  may  also  be  expressed  as  a percentage  error  PERRj  . 

PERR1  = 100  (e]  /DTOTj)  (256) 

If  PERR  is  sufficiently  small,  then  the  initial  allocation  by  steepest  ascent  is 

satisfactory  with  respect  to  total  damage  and  only  the  level  of  damage  to  the  defended 
target  subset  must  be  investigated.  If  PERRj  is  not  sufficiently  small,  then  a new 

allocation  must  be  made  using  the  procedure  outlined  below. 

The  Second  Allocation  - The  second  allocation  makes  use  of  the  same  X^’s  as  the 

initial  allocation.  The  general  procedure  is  to  find  the  target  type  for  which  there  is 
the  greatest  penalty  incurred  when  the  second-best  — rather  than  the  best  — weapon 
type  is  employed  against  it.  The  algorithm  to  determine  this  penalty  includes  the 
fraction  of  the  remaining  inventory  (initially  this  is  the  entire  inventory)  of  the  best 
weapon  (in  the  LaGrangian  profit  sense)  required  to  produce  this  profit  against  all  targets 
of  the  class,  FRACINV  , the  number  of  weapons  of  the  best  type  IWM  needed  for  a 


t 
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payoff  at  all,  ^+1 , the  best  LaGrangian  profit  per  target  for  target  type 

IT  , HMAXC^  and  the  second  best  LaGrangian  profit  for  target  type  IT,  HMAXN^  . 
Making  use  of  all  the  above,  the  penalty  parameter  for  target  type  IT,  CRIT-p  , is 
defined  as: 


CRITIT  = (NPj^  IT+  1)  • FRACINV(1.1  HMAXCIT  - HMAXN^)  (257) 


The  maximum  value  of  CRITjrp  over  the  target  set  is  then  selected  and  the  target 
type  which  produced  it  is  considered  the  most  critical  type  with  respect  to  the  allocation 

C RIT  = MAX  <C  RIT™  > (258) 

OVER  1 1 

ALL 

TARGET 

SUBCLASSES 


Provided  there  are  enough  weapons  of  the  best  type  for  use  against  at  least  a single 
target  in  the  most  critical  class  ITC  , these  are  applied  in  accordance  with  an  algorithm 
relating  (1)  the  total  shadow  value  of  the  NOTAS^  unassigned  weapons,  STOGO, 

(2)  the  total  shadow  value  of  all  weapons,  SVAL,  (3)  the  shadow  value  for  the  best 
weapon  X^^  » and  (4)  the  optimum  number  of  weapons  to  be  used  against  a single 

target  of  the  critical  type,  . The  number  applied  according  to  the 

algorithm  is: 


MAXWPS  = MAX  <[1  + (STOGO  - .3  SVAD/X^]  : WMM  ITC*> 

Where 

STOGO  = £ NOTAS^  XIW 

sval 


(259) 


(260) 


Now,  MAXWPS  of  the  type  IWM  are  applied  to  target  type  ITC  and  the  destruc- 
tion calculated.  NOTASjy^  and  STOGO  are  updated  to  reflect  the  expenditures.  In 

addition,  HMAXC^  and  HMAXN^  are  modified  as  necessary,  to  reflect  the  fact 

that  there  may  now  not  be  enough  weapons  if  type  IWM  was  the  best  or  second-best 
weapon  for  target  type  IT  . 

The  entire  procedure  is  now  repeated,  beginning  at  equation  (257)  until  STOGO  5 
,3  SVAL  . When  enough  weapons  have  been  allocated  to  critical  target  types  this  way  to 
make  STOGO  < .3  SVAL  , the  remainder  of  the  weapon  inventory  is  allocated  by  the 
steepest  ascent  method  used  for  the  initial  allocations. 


r ^ > #^SS«EZ3f7^ERr 


Testing  of  the  Second  Allocation  - The  second  allocation  is  testeu  using  the  total  expected 
value  destroyed  on  the  second  allocation,  DTOT2  , and  UBOUND  from  the  first  alloca- 
tion (Xp^’s  and  unchanged): 

e2  = UBOUND  - DTOT2 

PERR2  = 100  (e2/DTOT2)  (261) 

If  PERR2  is  sufficiently  small,  then  the  second  allocation  is  satisfactory  with 

respect  to  total  damage  and  only  the  level  of  damage  to  the  defended  target  subset  must 
be  investigated.  If  PERR2  is  not  sufficiently  small,  then  the  shadow  values  (Xp^’s) 

are  adjusted  using  the  procedure  outlined  below . 

Weapon  Shadow  Value  Adjustment  - If  the  second  allocation  outlined  is  not  satisfactorily 
close  to’ optimum,  the  optimum  criterion  is  adjusted  by  adjusting  the  weapon  shadow 
prices  (Xpy's) . This  is  carried  out  by  comparing  the  desired  inventory  of  each  weapon 

type  (as  derived  by  a heuristic  algorithm)  with  the  actual  inventory  of  the  weapon  type. 

If  more  weapons  of  type  IW  are  desired  than  currently  are  in  inventory,  the  price  is 
too  low  and  X^  should  be  increased.  Conversely,  if  fewer  of  weapon  type  IW  are 

desired  than  there  are  in  inventory,  X ^ should  be  decreased.  In  addition  to  changing 

the  individual  Xj^'s  , the  total  shadow  value  of  the  desired  inventory  should  match  the 

total  shadow  value  of  the  actual  inventory.  If  the  total  desired  inventory  shadow  value  is 
greater  than  the  actual  inventory  shadow  value,  it  indicates  that  the  imputed  weapon 
values  (X^,'s)  are  too  liigh  and  conversely.  When  the  desired  and  actual  inventories 

match  or  the  Xp^'s  are  sufficiently  close  and  the  desired  inventory  shadow  value  is 

sufficiently  close  to  the  actual  inventory,  the  Xpy's  are  used  to  re-evaluate  the  second 

allocation. 

Now  that  an  overview  of  the  procedure  has  been  given,  the  actual  steps  and  algo- 
rithms will  be  presented. 

The  first  step  is  to  determine  the  desired  inventory  for  each  weapon  based  on  the 
Xp^'s  derived  in  evaluating  the  second  allocation. 


The  inventory  increment  for  weapon  type  IW  against  target  type  IT  is: 

Axnimn,  _ WIW , IT^  MAX  ^ HIT,  IW  ' #95HIT^ 
ANWPDiw^it-  Hpp  ~ - 795Kit*> 


(262) 
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The  total  desired  inventory  NWPDRy  for  weapon  type  IW  is  then: 

NWPDIW-=  £ ANWPDIW  > IT  (263] 

Now  if  the  desired  inventory  for  each  weapon  is  compared  with  the  actual  inventory 
and  if  the  desired  inventory  is  less  than  the  actual  inventory  for  weapon  type  IW  , Xj^ 

is  reduced  by  5 percent. 


nwpdm<nwpiw 

A new  desired  inventory  is  now  computed  using  the  new  value  of  X ^ . Following 
this,  the  values  of  increments  to  the  X^' s , (AX^'s),  are  generated  using  a con- 
verging iterative  process.  Initial  values  for  the  AX^’s  are: 


^IW 

NWPDJW  > NWP. 

A ^iw  = r 

NWPDIW  £ NWP. 

Further  adjustment  is  then  accomplished  as  follows: 


1.3AX, 


NWPDIW<NWPIW 

and 

AXM  >0 

or 

NWPDjW  = NWPjyy 

or 

NWPDJW  > NWPrv 

and 

£> 

A 

O 

nwpdiw<nwprv 

and 

A ^ iyv  £ ® 

or 

NWPDjW  > NWP^ 

and 

AXm  so 

Further  refinements  are  now  made  to  the  AX^y's  based  in  part  on  the  absolute 


value  of  the  maximum  adjustment  made  so  far: 
XMAXL  = MAX^AXjyy  |> 


OVER  ' 
ALL 

WEAPON 

TYPES 


"v;_»y  #rvf  ■*>¥  *-  5vjj.  ' *****^^*«^P^ 


■S 


s 

\ 
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When  XMAXL  «J  .001  , (which  it  will  not  be  the  first  time)  sufficient  adjustment 
has  been  made.  If  not,  the  AX^'s  are  checked  and  adjusted  if  necessary  to  keep 

them  within  preset  bounds  in  two  steps  as  follows: 


Step  I Adjustment 


AXjw  - .8  , 


A X iw  < ” • 8 


AX^  -»XMAXL/20,  OsAX^  < XMAXL/20 


(268) 


Step  II  Adjustment 


IW 


A \m  -*  - XMAXL/20,  - XMAXL/20  < A X^  < 0 


AX^  -♦  .5 


A Xiw  > *5 


The  values  of  the  X ^'s  are  then  recalculated  taking  into  account  the  above 
AXjy^’s  and  a simultaneous  correction  factor  AX'  . Initially  AX’  = - .03  . 

^IW  **  ^IW^  + A X j\y)(l  A X ) 


(269) 


The  shadow  values  for  both  the  aggregate  desired  inventory  and  the  aggregate  actual 
inventory  are  now  computed  for  the  purpose  of  adjusting  AX'. 


(270) 


sval-Dx^nwp^ 

SVDL  ^XjwNWPD^ 

The  adjustments  to  AX'  are  now  made  similar  to  the  way  adjustments  were  made  to 
the  AX^y's  in  equation  (266): 


(271) 


AX'  -♦ 

- AX'/4 

SVDL  < SVAL  and  AX'  a 0 

SVDL  > SVAL  and  A X ' s 0 

t 

<3 

1.3  AX' 

SVDL  < SVAL  and  AX’  <0 

SVDL  > SVAL  and  A X ' > 0 

As  in  the  case  of  the  AX^'s  the  value  of  AX'  is  limited  as  follows: 
Step  I Adjustment 

AX’  - .8  AX’  < - .8 

AX'-*  XMAXL/20,  0 <;  AX'  < XMAXL/20 


(272) 


■»u  .-IV'  - 
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Step  II  Adjustment 

AX'*-  XMAXL/20, 
AX'  -»  .5, 


- XMAXL/20  < AX’  <0 
AX’  > .5 


The  above  sequence  is  now  repeated  beginning  at  equation  (262)  and  continued  until 
either  a present  number  of  iteration  (generally  100)  has  been  completed  of  XMAXL  < .001 
as  determined  in  equation  (267). 

Retesting  the  Second  Allocation  Based  on  Revised  Shadow  Values  - The  second  allocation 
is  now  tested  again  using  the  new  snadow  values  (Xp^vs)  just  derived.  The  procedure 

is  identical  to  the  procedure  followed  in  evaluating  the  initial  allocation  in  equations 
(249)-(256).  The  results  based  on  a new  value  of  UBOUND  are: 


e3  = UBOUND 


DTOTo 


(273) 

PERR3  = 100  (e3/DTOT2) 

Lf  PERR3  is  small  enough,  the  second  allocation  is  satisfactory  with  respect  to  total 

damage  and  only  the  level  of  damage  to  the  defended  target  subset  must  be  investigated. 

Lf  PERR3  is  not  small  enough,  then  a third  allocation  must  be  made,  by  means  of  the 

procedure  outlined  below. 

The  Third  Allocation  - The  third  allocation  is  generated  using  exactly  the  same  pro- 
cedure which  was  used  in  generating  the  second  allocation  but  utilizing  the  adjusted 
Xp^'s  discussed  earlier. 

Testing  the  Third  Allocation  - The  tiiird  allocation  is  tested  in  the  same  way  as  the 
second  allocation,  except  that  the  restriction  on  PERR^  for  a satisfactory  solution  is 

changed  from  1 percent  to  3 percent  (or  6 percent  if  it  is  not  desired  to  take  the  time 
necessary  for  3 percent  convergence).  If  PERR^  meets  the  criterion  with  respect  to 

total  damage,  then  only  the  level  of  damage  to  the  defended  target  subset  must  be 
investigated.  If  PERR^  does  not  meet  the  relaxed  criterion,  the  above  procedure  must 

be  repeated,  as  outlined  below. 

Additional  Adjustment  to  Shadow  Values  Re-evaluation,  Additional  Allocations  and  Tests 
ol  These  Allocations  - When  the  third  allocation  has  not  produced  satisfactory  results, 
the  Xp^'s  are  again  modified,  the  third  allocation  is  then  retested  and,  if  necessary, 

a fourth  allocation  is  generated  and  tested.  Similarly,  if  this  is  not  satisfactory,  the 
process  of  Xp^  adjustment,  re-evaluation,  new  allocations,  and  allocation  testing 

is  repeated  until  a satisfactory  allocation  is  found. 

Satisfying  the  Assured  Destruction  Constant  - When  a satisfactory  allocation  lias  been 
found  --  whe^ier  it  be  on  the  first,  second,  third,  etc,  allocation  — an  additional 
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constraint  must  also  be  satisfied.  That  constraint  is  the  need  to  destroy  at  least  some 
stipulated  fraction  (FRACP)  of  the  total  value  of  the  subset  of  defended  targets.  The 
fraction  actually  destroyed  is  termed  FRACPA.  This  is  also  known  as  the  preferred- 
target -destruction  criterion.  Thus,  starting  with  the  total  defended  target  value,  VPREF, 
and  the  total  defended  target  value  destroyed,  DPREF: 

VPREF  =£VIT 

SUBSET  OF 
DEFENDED 
TARGETS 


-.-DPREF  = £DIT 

SUBSET  OF 
DEFENDED 

targets 


(274) 


Then: 


FRACPA  = DPREF/VPREF 


(275) 


Now,  if  FRACPA  < FRACP  , all  non-defended  target  types  are  being  valued  too 
highly  and  a new  allocation  is  made . In  this  allocation  all  undefended  target  values  are 
multiplied  by  « = .003  to  reduce  their  values  on  this  next  allocation  to  .3  percent  of 

their  original  values . 

V.  .003  VTT 

1 1 target  1 ' (276) 

TYPE  IT 
UNDEFENDED 


The  allocation  procedure  is  the  same  as  that  previously  outlined,  beginning  with  a 
steepest  ascent  initial  allocation  and  proceeding  to  second,  third,  etc.,  allocations  as 
necessary.  The  only  difference  occurs  in  the  penalty  algorithm,  equation  (257),  in 
which  a multiplicative  factor  is  added.  The  value  of  this  factor  depends  on  whether  the 
allocation  just  previous  destroyed  enough  preferred  target  value. 

CRITFAC  = . 1 FRACPA  > FRACP 

CRITFAC  - 10  FRACPA  <;  FRACP  (277) 

CRITit  = CRITFAC(NPIWM  it+  1)  • FRACINV(1. 1 HMAXCIT  - HMAXNIT) 

After  the  allocation  procedure  is  again  completed  the  defended  target  value  destroyed 
is  recomputed  using  equation  (274)  and  the  fraction  of  total  defended  value  destroyed  is 
recomputed  using  equation  (275).  If  FRACPA  and  FRACP  are  now  within  0.02  of 
each  other,  or  if  FRACPA  < FRACP  (assured  destruction  is  inachievable),  the  allocation 
process  is  complete.  If  not,  the  undefended  target  value  multiplier  for  the  Nth 
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assured  destruction  iteration  is  varied  using  a linear  interpolation  method.  The  essen- 
tials of  the  process  are  explained  below  using  figure  18. 


FIG.  18:  RELATIONSHIP  BETWEEN  FRACTION  OF  DEFENDED  VALUE  DESTROYED 
AND  UNDEFENDED  TARGET  MULTIPLIER 

In  figure  18,  is  defined  as  the  highest  value  of  that  has  thus  far  produced  an 

acceptable  level  of  assured  destruction,  <xy  is  the  lowest  value  of  that  has  thus  far 

produced  an  unacceptable  level  of  assured  destruction.  Initially,  since  the  first  accept- 
able allocation  assumed  no  preference  for  defended -target  values,  «y  = 1 . Similarly 

since  the  first  iteration  using  = 0.003  was  successful  in  producing  an  acceptable 

level  of  assured  destruction,  then  initially  a = 0.003  . FRL  and  FRU,  then,  are 

L 

defined  as  the  fractions  of  defended  target  value  destroyed  (assumed  destruction  level), 
when  the  undefended  target  multiplier  takes  on  values  of  cc^  and  , respectively. 


Following  the  first  and  higher  iterations,  a^or  and  FRU  or  FRL  are 
updated  reflecting  the  value  of  FRACPA  obtained  on  the  Nth  iteration. 


*U  **  “N 
FRU  -»  FRACPA 


L N 

FRL  -»  FRACPA 


N a 1| 
N s 1 | 


f FRACPA  <FRACP 

I | FRACPA  - FRACPj  >.02 

r FRACPA  > FRACP 

I I FRACPA  - FRACP  | > 102 


(278) 
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Based  on  figure  18  and  the  above,  the  linear  interpolation  value  of  a ^ ^ for 
iteration  N+l  is  then: 

“N+i  = “l  + (aU  " “L)(FRL  ' pRACP)/(FRL  " FRU)  (279) 

The  value  of  undefended  targets  for  assured  destruction  iteration  N+l  are  then 
updated: 


VIT*  aN+lVIT  (280) 

TARGET 
TYPE  IT 
UNDEFENDED 

A new  allocation  is  then  processed  based  on  the  new  values  of  undefended  targets. 
When  this  allocation  is  completed,  it  is  tested  to  see  if  the  defended  target  fraction 
F RAC  PA  is  within  0.02  of  the  desired  fraction,  FRACP  . If  so,  the  allocation  is 
complete.  If  not,  but  the  limit  on  the  number  of  assured  destruction  iterations  has 
been  reached,  then  the  assured  destruction  iteration  allocation  in  which  FRACP  was 
exceeded  by  the  smallest  amount  is  taken  to  be  the  final  allocation.  If  FRACPA  is  not 
sufficiently  close  to  FRACP  (differing  by  more  than  0.02)  and  the  limit  for  assured 
destruction  iterations  has  not  yet  been  reached,  then  another  iteration  is  made  beginning 
with  the  updating  of  the  various  a's  in  equation  (278). 
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SECTION  III 


SUMMARY  OF  SUBROUTINE  OPERATIONS 
GENERAL  OVERVIEW 

As  a point  of  reference,  it  is  first  useful  to  review  the  typical  Code  50  scenario. 

While  many  other  possibilities  exist,  this  is  the  most  extensively  used.  fFirst,  it  is 
assumed  that  the  Soviet  Union  attacks  the  United  States  in  a counter -force -only  strike 
hitting  U.S.  land-based  missile  sites  and  non-alert  bombers  and  ballistic  missile  sub- 
marines. The  Soviets  are  assumed  to  hold  back  all  bombers  and  a smallTiumber  of 
missiles  for  the  third  strike,.  The' Soviet  attack  is  optimized  to  maximize  the  U.S. 
weapon  value  destroyedr  After  the  destruction  is  taken  into  account,  all  surviving  U.S. 
forces  attack  the  Soviet  cities  in  an  attack  that  is  optimized  to  maximize  fatalities, 
subject  to  the  assured  destruction  constraint  that  X percent  of  the  value  in  terminally 
defended  cities  must  be  destroyed.  After  the  second  strike,  the  Soviets  make  a third 
strike  against  U.S.  cities  with  their  remaining  weapons,  to  kill  a maximum  number  of 
people  ft  On  this  basis,  the  discussion  can  now  proceed  into  how  Code  50  operates  to 
work  tTTrough  this  scenario. 

The  Code  50  program  is  divided  into  two  parts.  The  first  consists  of  the  war-game 
model  itself  and  the  data  input  systems.  This  segment  of  Code  50  reads  in  data,  initial- 
izes for  the  ensuing  computations,  calculates  single-shot  kill  probabilities,  analyzes 
weapon -interceptor  engagements  for  both  missiles  and  bombers,  and  in  general,  prepares 
the  weapon  and  target  data  for  the  second  part  of  the  program  which  will  ietermine  the 
weapon-to-carget  allocation.  The  various  components  of  this  part  of  Code  50  are  Code 
50,  the  main  executive  program;  GENX  and  DATAX,  the  input  routines;  STRIKE,  a 
secondary  executive  and  general  calculating  routine;  MBPEN,  a single-island  random 
area  penetration  subroutine  for  both  missiles  and  bombers;  SSKCALC,  which  calculates 
single-shot  kill  probabilities;  ant!  MINNUM,  which  computes  the  minimum  number  of 
weapons  of  each  type  which  must  be  invested  at  each  type  of  target  in  order  to  obtain  a 
payoff. 

The  game  model  segment  of  Code  50  is  interfaced  with  the  weapon-to -target  allocation 
through  subroutine  TRANS.  Control  then  passes  to  subroutine  ASSIGN  which  is  the  execu- 
tive subroutine  for  the  allocator  section  of  Code  50.  When  the  allocation  is  complete, 
the  damage  is  computed,  control  is  returned  to  Code  50,  the  attacker  becomes  the  de- 
fender and  vice-versa,  and  subroutine  STRIKE  is  called  again.  The  procedure  is  carried 
out  three  times  or  less,  although  the  initializing  and  data -reading  subroutines  are  only 
called  on  the  first  strike. 

THE  WAR  GAME  SEGMENT 

The  basic  purpose  of  the  war  game  segment  of  Code  50  is  to  feed  to  the  inside  (or 
allocator)  segment  lists  of  available  weapons  to  allocate  to  the  various  targets  along  with 
all  requisite  information  concerning  both  weapons  and  targets.  The  subroutine  calling- 
sequence  for  the  war  game  is  as  shown  in  figure  19,  the  events  represented  by  this 
sequence  are  as  follows.  The  executive  routine  Code  50  sets  the  attacker  and  defender 
indices  and  calls  the  data  input  subroutines,  GENX  and  DATAX.  These  read  in  informa- 
tion and  also  compute  some  parameters  such  as  the  number  of  weapon  and  target  types 
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for  each  side.  Control  then  returns  to  Code  50  where  an  initializer  transfers  data  from 
the  fixed  arrays  to  the  working  arrays.  Then  a secondary  executive  routine,  STRIKE, 
is  called.  It  sets  the  number  of  weapon  classes  and  target  classes  for  attacker  and 
defender  for  later  use.  It  then  obtains  the  target  values.  There  is  an  option  to  either 
read  in  target  values  or  calculate  them.  The  calculation,  of  course,  applies  only  to 
weapons  considered  as  targets  since  the  value  of  a city  has  been  taken  as  its  population. 
The  other  task  of  STRIKE  is  to  set  flags  in  arrays  to  indicate  missiles,  targets  with 
extended  areas  to  which  the  square -root  law  is  applied  to  calculate  survival,  targets 
which  are  terminally  defended,  and  to  update  the  weapons  left. 

The  subroutine  that  calculates  the  area  penetration  probabilities,  MBPEN  is  now 
called.  The  bomber  and  missile  area  penetration  models  covered  in  section  II  are 
included  in  MBPEN.  MBPEN  is  called  with  an  argument  so  it  can  be  used  for  either  side 
at  any  time,  not  merely  the  attacker.  This  capability  is  required  if  one  assumes  that 
the  value  of  a weapon  as  a target  is  in  part  a function  of  its  delivered  yield.  Thus,  given 
the  index  for  one  side,  the  first  step  is  to  find  the  other  index.  MBPEN  then  proceeds  to 
find  the  number  of  bombers  of  each  type  arriving  at  the  area  defenses  as  well  as  the  total 
number.  There  is  an  option  to  read  the  kill  probability  of  a bomber  by  a fighter  instead 
of  calculating  it.  The  penetration  probability  is  defined  as  the  probability  that  a bomber 
is  not  killed  before  delivery  multiplied  by  its  ASM  inflight  launch  reliability.  This  last 
is  unity  for  gravity  bombers.  The  last  point  here  is  the  calculation  of  NWP  which  is  the 
number  of  reliable  independently  targetable  weapons. 

Subroutine  SSKALC  calculates  the  single -shot  kill  probabilities  and  modifies  them 
to  include  any  restrictions  placed  on  the  weapons  such  as  not  sending  bombers  against 
alert  missiles.  This  is  done  using  the  single -shot  kill  models  discussed  in  section  II. 

Subroutine  MINNUM  gives  the  minimum  number  of  weapons  of  a given  type  required 
to  obtain  a payoff  at  a target.  The  bomber  and  missile  terminal  penetration  models  dis- 
cussed in  section  II  are  used  for  this  purpose. 

This  completes  the  description  of  the  game  model  and  how  the  specific  game  models 
of  section  II  fit  together.  To  recapitulate,  this  section  of  the  program  feeds  to  the  inside 
a list  of  numbers  of  reliable  weapons,  kill  probabilities,  number  of  weapons  required  for 
a payoff,  a target  list,  and  numbers  of  weapon  and  target  types.  At  this  point  the  weapon - 
to-target  allocation  routine  takes  over.  Then  after  the  allocation  is  complete,  and  a 
summary  result  is  printed,  control  returns  to  Code  50  where  the  attacking  and  defending 
sides  are  exchanged,  the  strike  number  is  increased  by  unity,  and  the  sequence  starts 
again. 

THE  WEAPON  ALLOCATION  SEGMENT 

The  calling  sequence  for  the  weapon  allocation  segment  is  handled  by  subroutine 
ASSIGN.  Its  flow  chart  is  shown  in  figures  20-23.  First,  it  calculates  several 
quantities  including  the  total  and  terminally  defended  target  system  values.  Next  it  calls 
an  initializing  routine  (INITIAL),  and  then  performs  the  steepest  ascent  initial  laydown 
allocation.  As  discussed  in  section  II  in  the  Weapon  Target  Allocation  Model,  this  alloca- 
tion is  a compormise  in  that  it  neither  lays  down  all  of  a given  weapon  type  at  once,  nor 
does  it  lay  the  weapons  down  individually.  Rather  it  makes  the  laydown  using  one-fifth 
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of  each  weapon  type  on  the  best  targets.  The  best  target  for  a given  weapon  type  is  found 
in  subroutine  MAXPAY.  As  some  number  of  targets  of  a given  type  have  weapons  assigned 
to  them,  new  target  classes  are  created  for  the  targets  of  that  type  which  have  not  as  yet 
been  touched.  This  is  done  using  subroutine  SPLIT.  When  the  laydown  is  completed,  one 
more  weapon  of  each  type  is  then  assigned  to  the  best  target  for  that  weapon.  The  bare 
payoff  or  best  alternative  use  for  a single  weapon  is  taken  to  be  the  starting  value  of  the 
LaGrange  multiplier.  Having  obtained  a set  of  multipliers,  the  next  step  is  to  evaluate 
the  allocation  just  made.  To  do  this,  one  must  have  in  some  sense  an  optimum  allocation. 
This  optimum  is  obtained  using  subroutines  MAXLAG  and  EVAL.  The  first  of  these  finds 
the  maximum  LaGrangian  payoff  for  each  target  if  the  optimum  number  of  the  best  weapon 
type  is  used  at  each  target.  Subroutine  EVAL  then  compares  the  total  payoff  using  these 
optimal  LaGrangian  payoffs  with  the  actual  payoff  obtained  using  the  laydown  allocation. 
Generally,  the  fractional  error  will  be  over  the  preset  limit  so  a new  allocation  must  be 
made.  This  new  allocation  is  made  on  the  most  critical  targets  first.  These  are  the 
targets  on  which  the  largest  differences  are  obtained  in  the  LaGrangian  payoff  if  the 
optimum  weapon  were  out  of  stock.  This  allocation  proceeds  until  seventy  percent  of  tne 
total  snadow  value  of  tne  weapons  has  been  expended.  The  remaining  thirty  percent  is 
allocated  using  the  steepest  ascent  laydown  algorithm.  Again  the  allocation  is  evaluated. 

If  it  is  not  good  enough,  the  LaGrange  multipliers  aie  adjusted  using  subroutine  ADJLAM. 
This  subroutine  first  calls  OPLAG  which  finds  for  the  given  types  of  weapons  in  stock 
the  optimum  number  of  each  weapon  to  have.  The  difference  between  the  optimum  num- 
ber and  the  actual  number  is  the  criterion  by  which  the  lambdas  are  adjusted.  This 
adjusting  involves  a number  of  heuristic  considerations  and  is  both  complicated  and  time- 
consuming.  After  it  is  completed,  another  allocation  is  made.  Experience  has  indicated 
that  this  allocation  is  almost  always  very  close  to  the  optimum,  although  additional 
allocations  may  need  to  be  made  in  some  cases  using  the  loop  mechanism  shown  in 
figures  20-22.  This  completes  the  weapon -to -target  allocation.  If  there  is  no  assured 
destruction  required,  ASSIGN  calls  CONSOLID  to  find  any  repetitions  in  the  allocation 
list  and  to  combine  such  repetitions  into  a single  target  type  subclass.  ASSIGN  then 
makes  a final  call  on  EVAL,  calculates  damage  and  calls  on  TARGSUM  to  print  out  the 
damage.  Then  ASSIGN  calls  WPSUM  to  print  out  a summary  of  the  attacking  side  weapons. 
Finally,  PRALLOC  is  called  to  print  out  a summary  of  the  allocation.  Control  then 
returns  to  Code  50. 

On  the  other  hand,  if  there  is  some  assured  destruction  requirement,  a somewhat 
longer  calculation  is  made.  A ciieck  is  first  made  to  see  if  the  assured  destruction  was 
obtained  on  the  optimum  allocation.  If  so,  the  allocation  is  complete.  If  not,  the  value 
of  all  the  non-preferred  targets  is  set  to  zero,  and  we  start  over.  On  returning  to  this 
point,  a check  is  made  to  see  if  it  is  ever  possible  to  obtain  the  assured  destruction  on 
the  preferred  targets.  If  not,  a message  is  printed  and  the  allocation  is  complete.  If 
the  assured  damage  level  was  obtained  with  the  non-preferred  target  values  set  to  0.3 
percent  of  their  values,  an  iterative  technique  is  used  as  explained  in  section  II  to  find 
the  factor  by  which  the  non -preferred  target  values  are  to  be  multiplied  to  obtain  as 
nearly  as  possible  the  exact  assured  damage  level.  The  allocation  is  then  returned  to 
Code  50  through  MINNUN,  TRANS  and  STRIKE,  as  shown  in  figure  19. 
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FIG.  19:  OUTLINE  OF  CNA  50  OUTER  CALCULATION 
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FIG.  20:  OUTLINE  OF  INNER  CALCULATION: 
FIRST  ALLOCATION  CALCULATION 
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FIG.  21:  OUTLINE  OF  INNER  CALCULATION: 
LATER  ALLOCATION  CALCULATION 


I 

I 


RECOMPUTE 

Xiw's 

(MAXPAY) 


i 


I 


\ 

\ 

\ 


FIG.  22:  OUTLINE  OF  INNER  CALCULATION 
ADJUSTING  THE  LAGRANGE  MULTIPLIERS 


FIG.  23:  OUTLINE  OF  INNER  CALCULATION: 
REQUIRED  DESTRUCTION  CALCULATIONS 
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